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Abstract 



Different aspects of the non-equilibrium chiral phase transition as well as in- 
medium properties of the scalar and pseudoscalar mesons is investigated using 
quantum field theoretical descriptions at nonzero temperatures and densities. The 
chiral phase transition at nonzero temperature and baryon chemical potential is 
studied at mean field level in the sigma model that includes quark degrees of 
freedom explicitly. For small bare quark masses the critical point separating the 
first order phase transition line and the smooth crossover region is determined. 
The behavior of quark and meson masses with changing temperature and chemical 
potential is examined. Adiabatic lines are computed showing that the critical point 
does not serve as focusing point in the adiabatic expansion. The linear sigma model 
without quarks but with field fluctuations beyond the mean field is used to derive 
coarse-grained equations of motion for the soft homogeneous and inhomogeneous 
chiral condensate fields coupled to a heath bath. Multiple effects of the heat bath 
on the chiral condensate fields is studied based on linear response theory. The 
self-consistently solved gap equation is used to determine the order of the phase 
transition and the temperature-dependence of meson masses. Using these masses 
the decay widths of pions and sigma mesons are computed at nonzero temperature 
at one- and two-loop order in perturbation theory. Numerical results show that in 
the phase transition region not only is the damping of sigma mesons significant, 
but also that of the pions. The in-medium modification of the pion dispersion 
relation is examined. A deviation at finite temperature of the velocity of massless 
pions from the speed of light is found. Separate investigation is carried out for the 
model with exact and explicitly broken chiral symmetry, making sure Goldstone's 
Theorem is fulfilled when chiral symmetry is spontaneously broken. 
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Chapter 1 
Introduction 



According to the Big Bang Theory, the Universe started from an extremely hot and 
dense singularity. As it expanded, it cooled rapidly, passing through several phase 
transitions. The GUT (Grand Unified Theory) transition at about 1CT 34 seconds 
after the Big Bang, and at a temperature of about 10 29 K was followed by the 
EW (Electro- Weak) transition at about lCT 12 s and 10 15 K. When the Universe was 
about 1CT 6 seconds old, and the temperature was about 150 — 200 MeV ~ 10 12 K 
the confinement phase transition took place. This has great significance since it 
is due to this transition that the world we live in, the hadronic world, has been 
formed. 

The theory that governs the nature of the confinement transition is QCD 
(Quantum Chromodynamics), the fundamental theory of strong interactions. QCD 
is a relativistic quantum field theory and as such very complicated. Its finite tem- 
perature properties cannot be determined analytically. Therefore, instead of direct 
calculations, physicists need to rely on other approaches. Two major trends are 
lattice QCD and phenomenological effective field theories. Throughout this disser- 
tation the focus is on effective field theories. 

Why is it so difficult to analyze QCD directly? The usual approach when 
studying a field theory is to perform a perturbative expansion in the coupling, 
keep the dominant terms of the expansion, and then look at the ground state 
and the spectrum of excitations. Such perturbative treatment is reasonable in the 
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asymptotically free region of QCD where the coupling is weak, a s <l, for large 
momentum transfers (q 2 ^> (lOOMeV) 2 ) or small distances 0, 0. In this region 
quarks and gluons, the basic degrees of freedom in QCD, interact weakly. In our 
world, though, quarks and gluons appear only in confined states. Confined by the 
strong interaction, they form hadrons (mesons and baryons). At hadronic energy 
scales (q 2 < (lOOMeV) 2 ), or large distances, the QCD coupling constant is large, 
a s ^> 1, making a perturbative expansion not feasible. 

As an important feature of QCD, it has been predicted and now widely ac- 
cepted that at temperatures about 150 MeV a deconfinement of hadrons happens 
according to which quarks and gluons move freely. Here perturbative approaches 
should be quite reasonable. When studying QCD at finite temperature, however, 
one encounters infrared divergences which causes the breakdown of perturbation 
theory ||. A rather successful non-perturbative method is putting QCD on a 
lattice and performing numerical analysis. The studies done in the pure gauge 
sector (only gluon degrees of freedom) are of great importance for understanding 
high temperature QCD. Developments including fermion dynamics (quarks) have 
begun, but improvement of the present algorithms is required. The most recent pa- 
pers by the lattice-community H || report promising results at finite temperature. 
An alternative way to study the fundamental theory is to construct effective field 
theories. These incorporate some aspects of QCD and ignore others. Therefore, 
by sacrificing some of the richness of the full theory, one gains some simplifications 
that in certain cases allow even for analytic solutions. Examples of such theories 
are the SU(2) R x SU(2) L and SU(3) R x SU(3) L § linear sigma models, the 
nonlinear sigma model ||12||, and the Nambu-Jona-Lasinio model H. 



Another important characteristic of the QCD structure is the non-vanishing 
expectation value of quark operators in vacuum. Present at low temperatures, 
the non-zero chiral condensate is nothing else but the Bose condensate of the 
underlying quark degrees of freedom, much like Cooper-pairs in superconductors, 
and liquid 3 He and electrons in superfluids [|l0j- In condensed matter physics the 
normal state does not include condensates. These can develop after crossing a 
phase transition at a characteristic temperature T c . In nuclear physics this is the 
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opposite, since the condensate is present in the normal state. It is widely believed 
that above a critical temperature of T c ~ 150 MeV it is possible to generate regions 
where the condensate is dissolved. 

What is the underlying physics that governs the formation and disappearance 
of condensates? The answer, common in Nature and well-rooted in most major 
parts of physics, is symmetry and breaking of symmetry. These concepts, first 
introduced for explaining superconductivity and superfluidity 0], are now widely 
applied in particle and nuclear physics. It is common to all of these situations 
that a continuous symmetry is spontaneously broken at low temperatures (ordered 
phase) and subsequently restored at high temperatures (disordered phase) by go- 
ing through a phase transition. Lattice results suggest [|13j that the transition 
between symmetric and broken phases happens together with the hadronization, 
the confinement of quarks and gluons into hadrons. 

A phase transition that happens in thermal equilibrium can always be char- 
acterized by an order parameter which is nonzero in one phase and vanishes in 
the other. For details the reader is referred to QCD with two massless 

quarks possesses chiral symmetry, described by the SU(2) R x SU{2) L group. This 
is isomorphic to the 0(4) group which describes the rotational symmetry in fer- 
romagnets. In both cases the continuous symmetry is spontaneously broken at 
low temperatures, resulting in the emergence of Goldstone bosons [Q. These are 
spin- waves in ferromagnets and massless pions in the two flavor QCD. The order 
parameters are the magnetization and the chiral condensate, respectively. If the 
transition is of second order universality arguments become rather powerful, pro- 
viding qualitative features of the two massless flavor QCD phase transition ||36|| . 
Nature, however, supplies a different situation. First of all, quark masses are non- 
zero. The small quark masses explicitly break chiral symmetry, just as an external 
magnetic field breaks the rotational symmetry of a ferromagnet. Explicit sym- 
metry breaking then alters the nature of the phase transition. The second order 
transition in a ferromagnet becomes a smooth cross-over. The change in the order 
of the chiral phase transition is non-trivial (see Chapter ^| and related references) . 
Second, at temperatures of interest three flavor QCD studies become necessary. 
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The effect of strange quarks on the phase transition has been recently discussed 
in [15 1 and references therein. Third, the phase transition need not happen in lo- 



cal thermal equilibrium circumstances. In case of an out-of-equilibrium transition 
there are no universality arguments one can rely on. However, such transitions are 
known in condensed matter physics (for example in superfluids and liquid crys- 
tals) and can guide us in nuclear physics. There is a major difference, though: 
while classical approach is sufficient for describing a phase transition in condensed 
matter physics, quantum field theoretical treatment is required in nuclear physics. 
For recent advances in equilibrium and non-equilibrium aspects of quantum field 
theories at finite temperature the reader is referred to |T7[ . 

In this dissertation different aspects of the non-equilibrium chiral/confinement 
phase transition as well as in-medium properties of the scalar and pseudoscalar 
mesons are investigated within an effective quantum field theoretical model, the 
linear sigma model. This thesis is structured as follows: In Chapter |2| the chiral 
phase transition at nonzero temperature and baryon chemical potential is stud- 
ied within the linear sigma model including quark degrees of freedom explicitly. 
In Chapter ||] we derive equations of motion for the long-wavelength chiral con- 
densate fields which are in contact with a heat bath. Chapter |] presents the 
temperature-dependence of the equilibrium value of the condensate and that of the 
self-consistently calculated meson masses. In Chapter [| the effect of the thermal 
medium on the velocity of long-wavelength pion modes is discussed and different 
dispersion relations are examined. The topic of Chapters || and [7] is the dissipation 
of long-wavelength sigma and pion condensate fields at one- and two-loop order, 
respectively. 



1.1 Chiral Symmetry in Quantum Chromody- 
namics 

A symmetry in Nature arises whenever a change in a variable leaves the essential 
physics unchanged. This translates to field theory as invariance of the Lagrangian 
under field transformations. Chiral symmetry is characteristic of theories with 
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massless fermions, and as such it is a fundamental symmetry of QCD with two 
massless quarks. Of course, QCD describes the physics of all six quark flavours. 
Nevertheless, the low energy limit of QCD is dominated by u and d quarks. The 
discussion presented below can be naturally generalized from two flavors to three, 
provided the s quark is considered to be light as well (although this can easily be 
seen as unreasonable at very low temperatures). Other flavors do not play a role 
below the natural scale of QCD, Aqcd ~ 100 MeV, at which the theory becomes 
strongly coupled. Mostly in regimes with temperatures T 3> Aqcd and/or baryon 
density % > Aq CD , where QCD describes quarks and gluons that are weakly 
interacting, is where the physics of heavier quarks becomes relevant. 
The massless QCD Lagrangian is 

Lqcd = qirf U + W^t) 1 ~ \ F ^ U , (1.1) 

where q is an Nf- dimensional Dirac spinor describing the quark field, Nf being the 
number of flavors. The index a refers to the color of the gauge field A a , representing 
the gluons. The Yang-Mills field strength tensor is 

F; v = d,A a v - d v Al - 9t hc A\A c u . (1.2) 

Using the chirality projection operators 

qn.,L = \(l±l*)q (1-3) 

the quark fields can be decomposed into their right- and left-handed components. 
Then the fermionic part of Lagrangian (|1 . 1| ) , rewritten as 

Lqcd = <ZW (dp + ig y^j q R + qii<f (d, + ig^A^j q L , (1.4) 
is invariant under the global transformations 

q R ,L -> e- iaR ^q R ,L • (1.5) 



This basically corresponds to rotating independently the right- and left-handed 
components of the spinors that describe the quark fields and are called chiral trans- 
formations. The symmetry group of this transformation is SU(Nf)n x SU(Nf)L. 
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With every invariance of the Lagrangian under a global transformation there is 
a conserved current associated with it. This is Noether's Theorem. The Noether 
current corresponding to the chiral transformations is the axial vector current 

= QrI^Qr- QlI^Ql ■ (1-6) 

This current is not only conserved when the quarks are massless, 

= , (1.7) 

but also does not connect spinors with different chiralities. Dirac bilinears like 
qq and qy'q have the same quantum numbers as scalar and pseudoscalar mesons, 
respectively. The vacuum expectation value of the scalar bilinear, what we refer 
to as the scalar condensate, is nonzero 

(qq) = (q L q R ) + o . (1.8) 

This signals that the vacuum mixes left- and right-handed quarks, allowing quarks 
to acquire an effective mass as they move through the vacuum. We say that the 
symmetry is spontaneously broken. This means that the Lagrangian maintains 
invariance under the transformations but the ground state breaks that invariance. 
For two flavors we say that the global chiral symmetry SU (2) R x SU (2) l is sponta- 



neously broken to a subgroup SU(2) L+R . According to Goldstone's Theorem JT4 
the result of spontaneous symmetry breaking is the existence of massless bosons. 
For two massless quark flavors there exists three massless pseudoscalars. These 
are the three pions. 

The massless quark situation is an idealization. One has to account for small 
quark masses. The Lagrangian then contains a mass term of the form 

— qmq (1.9) 

where m is a mass-matrix. The presence of this term results in the mixing of left- 
and right-handed quarks under the chiral transformation (|1.5| ): 

qLmq R + q R mq L (1.10) 
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Thus chirality is not preserved anymore. Because the symmetry is not exact the 
axial current is not conserved either 



8^f b = 2imqj 5 q . (1.11) 

This non-vanishing divergence measures the degree to which chiral symmetry is 
broken. The small pion mass is the result of explicit symmetry breaking by small 
quark masses. A rather comprehensive discussion of this physics is presented in a 



recent article [21| 



1.2 The Linear Sigma Model 

Effective field theories are extremely useful tools in understanding QCD physics. 
In such theories the QCD Lagrangian ( II . 1| ) given in terms of the basic quark and 
gluon degrees of freedom, is exchanged for Lagrangians built out of confined states 
of these (baryons, mesons, glueballs). The soundness of this shift in the degrees of 
freedom is demonstrated in |Tj| using the path integral formalism. For details of 
this procedure the reader is referred to [|1^, 

The linear sigma model |7| is one of the most instructive of all effective field 
theories. It describes the low energy non-perturbative sector of QCD extremely 
well. When using this model one gains access to analytical derivations and acces- 
sible numerical solutions, albeit at the expense of losing some properties, such as 
confinement. The Lagrangian of the sigma model is 

C = i> [z 7 ^ - g(a + z 75 f ■ tt)] V + \ (d^a + <9 M vr • d^) - U(a, tt) . (1.12) 

Here ip is a fermion field which, in early works has been considered to describe 
nucleons, but can also be identified with quarks as we shall discuss this in Chapter 
0. The scalar field a and the pseudoscalar field 7? = (71 - !, 7r 2 , 7r 3 ) together form a 
chiral field $ = (a, 7?). The usual choice for the potential is 

U(a,n) = -(a 2 + ti 2 - v 2 ) 2 - Ha . (1.13) 
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The coupling constant between mesons is A and between mesons and fermions is g 



and are dimensionless. The theory thus is renormalizable. The Lagrangian ( 1.12 
is invariant under chiral 577(2) l x SU (2) r transformations if the explicit symmetry 
breaking term Ha is zero. The parameters of the Lagrangian are usually chosen 
such that chiral symmetry is spontaneously broken in the vacuum. The potential 
resembles the bottom of a wine bottle. The usual choice for the expectation values 
of the meson fields are (a) = v = f n and (tt) = 0, where f n = 93 MeV is the pion 
decay constant. Different choices of ground states can be interchanged by a chiral 
rotation in the circle that is the bottom of the potential. Defining fluctuations 
about the vacuum as 

a = v + a' (1.14) 



and inserting this into Q1.12 ) results in a mass term for fermions M = gf n and for 
the sigma m 2 = 2\f%. These mass terms break the symmetry of the Lagrangian, 
illustrating that spontaneous breaking of symmetry corresponds to a particular 
choice of ground state out of the degenerate set. Associated with this mechanism 
are massless pion excitations. 

The symmetry is explicitly broken for any H ^ 0. The vacuum is not degen- 
erate anymore, the potential is tilted toward a unique value, and pions are not 
massless. The constant H is fixed by the PCAC relation which gives H = f n m1, 
where m n = 138 MeV is the pion mass in vacuum. Then one finds v 2 = f^ — m^/X. 
The coupling constant A is determined by the sigma mass, m 2 = 2Xf 2 + m 2 T , which 
is about 600 MeV in vacuum, yielding A ~ 20. The coupling constant g is usually 
fixed by the requirement that the fermion mass is the mass of the nucleon or that 
of the constituent quark in vacuum, which is about 1/3 of the nucleon mass. 

The linear sigma model was studied quite intensely at finite temperatures (see 
the pioneering works & ffl). The isomorphism of the 577(2) x SU(2) symmetry 



group with the 0(4) group led to studies based on universality arguments |36 
Consequently, the order parameter of this model, the quark condensate (qq), and 
its temperature dependence is expected to have similar behavior as magnetization 
in a ferromagnet. The sigma model without including fermions undergoes a second 
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order phase transition which leads to the restoration of broken symmetry above a 
critical temperature [71 . 



1.3 Ultra- Relativist ic Heavy Ion Collisions 

Ultra-relativistic heavy ion collisions offer a possibility to study the properties and 
the behavior of high temperature hadronic matter with the hope of detection of the 
deconfmement phase transition. Experimental programs at the European Center 
for Nuclear Research (CERN) SPS and Brookhaven National Laboratory (BNL) 
AGS have focused on this research up until 1999. In 2000 the Relativistic Heavy 
Ion Collider (RHIC) at BNL started running at never before reached energies. 
There are plans for even higher energies at the LHC at CERN, which currently is 
under construction, and will starts its operation in 2007. 

In these collisions heavy atoms such as lead (Pb 208 ) and gold (Au 197 ) are 
stripped of their electrons, accelerated to relativistic speeds and smashed into each 
other. In central collisions of two fast ions a considerable compression of matter 
occurs, leading to a very hot and dense nuclear matter. After such energetic im- 
pact (100 GeV/nucleon in the center of momentum frame at RHIC) temperatures 
are expected to reach above the critical temperature (~ 150 MeV), creating a de- 
confined form of hadronic matter. In this state of matter quarks and gluons are 
moving freely without the constraint of a confining potential. This soup of the 
basic constituents of matter is known as quark-gluon plasma (QGP). The possible 
formation and evolution of QGP is the most stimulating phenomena in relativistic 
heavy ion physics and has been the subject of intense investigations in the last cou- 
ple of decades both from theoretical and experimental perspectives. For a review 



on the advances see [22 



Due to the very high internal pressure the system expands and cools rapidly. 
The initial stages of a QGP are dominated by hard scatterings of quarks and 



qluons, leading to thermalization |23[ . Once a local thermal equilibrium is achieved 
a hydrodynamical description is feasible [^3J . During this hydrodynamic evolution 
the plasma expands and cools adiabatically. When reaching a critical temperature 
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rehadronization sets in. After the newly formed hadrons (mostly pions) freeze-out 
(cease interacting), these fly to the detectors. 

There is a wast amount of physics originating in the above very schematically 
described scenario. This is partly because only indirectly, by inferring back from 
what is recorded in the detectors, can the formation of QGP be detected (see for 
§). For any further phenomeno logical or technical details, or for 



example [45 



viewing images of current heavy ion experiments at RHIC see [25 
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Chapter 2 

Chiral Phase Diagram 



In this chapter the order of the chiral symmetry restoring phase transition is stud- 
ied at finite temperatures and baryon densities. The linear sigma model coupled 
to quarks is used in the mean field approximation to verify the existence of the 
tricritical point and to find the line of first order transition in the phase diagram. 
Then the behavior of quark and meson masses with changing temperatures and 
chemical potentials is examined. Finally, the adiabatic lines are computed in the 
(T, hb) plane. 



2.1 Review of Previous Results 

The idea that the spontaneously broken chiral symmetry of the QCD vacuum is 
restored at high temperatures is well accepted. The order of the transition, though, 
has not been determined unambiguously. 

It is important to determine the order of the chiral transition, as this influences 
the dynamical evolution of the system. It has been shown that a first order tran- 
sition in rapidly expanding matter may manifest itself by strong non-statistical 
fluctuations due to droplet formation [p8| . In the case of strong supercooling it 



may lead to large fluctuations due to spinodal decomposition |29|, ^Of. In a sec- 
ond order phase transition one may expect the appearance of critical fluctuations 
due to a large correlation length f3l| . Experimentally, large-acceptance detectors 
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are now able to measure average as well as event-by-event observables, which in 
principle can distinguish between scenarios with a first order, a second order, or 
merely a crossover type of phase transition. 

Theoretically, the QCD phase diagram in the (T,/ib) plane has recently re- 
ceived much attention |3l], [33], |34|, [35|1 . QCD with two flavors of massless quarks 
has a global SU{2) L x SU{2) R symmetry (see Chapter [I]). This symmetry is 
spontaneously broken in the QCD vacuum, such that the order parameter ((IlQr) 
acquires a non-vanishing expectation value, where q is the two-flavor quark field. 
At zero baryon chemical potential the effective theory for this order parameter is 
the same as the 0(4) model which has a second order phase transition. Therefore, 
by universality arguments [36|, the chiral transition in two flavor QCD is likely to 
be of second order at = 0. Nonzero quark masses introduce a term in the QCD 
Lagrangian which explicitly breaks chiral symmetry. Then there is no theoretical 
argument for the existence of a phase transition in the strict sense. In this case 
the second order transition becomes a crossover. 

At nonzero baryon chemical potential it is more difficult to infer the order of 
the chiral transition from universality arguments |p7[| . One commonly resorts to 
phenomenological models to describe the chiral transition in this case. Depending 
on the parameters of these models they may predict a first order, a second order, 
or a crossover transition. However, if there is a s econd order phase transition for 
fiB = and nonzero T and a first order transition for small T and nonzero fis, 
then there exists a tricritical point in the (T, hb) plane where the line of first order 
phase transitions meets the line of second order phase transitions. For nonzero 
quark masses, this tricritical point becomes a critical point. Very recent lattice 
calculations |27| reported the determination of this critical point. 

It has been proposed |3l| that this point could lead to interesting signatures in 
heavy-ion collisions at intermediate energies if the evolution went through or close 
to this critical point. At this point, susceptibilities (e.g. the heat capacity) di- 
verge, the order parameter field becomes massless and consequently exhibits strong 
fluctuations. This could be detected in event-by-event observables. A thorough 



discussion on how to find the QCD critical point is given in |32 |. 
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2.2 Thermodynamics 



The thermodynamics of two effective field theoretical models, the linear sigma 
model and the Nambu-Jona-Lasinio model, has been recently discussed and com- 
pared in Here we present the details and results of the investigation of the 
linear sigma model which includes quark degrees of freedom. The use of such 
model is justified since high temperatures are dominated by quarks and low tem- 
peratures are dominated by mesonic excitations. Our model differs from other, 
more common realizations of the sigma model, where quark degrees of freedom 
are neglected but mesonic excitations are included (see for instance reference [p9|). 
We believe that our approach is more justified at high T and fi when constituent 
quarks become light but mesonic excitations are heavy (see below). Of course, 
due to the confining forces, at low T and \x quarks and antiquarks will recombine 
into mesons, baryons and antibaryons. We can only hope that this hadronization 
process will not change drastically the character of the chiral transition which we 
study here. Let us recall the Lagrangian ( |1.12| ) of this model 

C = q [i 7 "fy - g(a + i l5 r ■ 7?)] q + i (d.ad^a + drf ■ 0"tt) - U(a, 7?) , (2.1) 

where the potential is 

U(a,n) = ^ (a 2 + tt 2 - v 2 f -Ha . (2.2) 

Here q is the light quark field q = (u,d), the scalar a field describes the sigma 
meson and tt = (7Ti, iT2, 713) the three pions. The choice of the parameters of 
the Lagrangian are described in section L2. When H = chiral symmetry is 
spontaneously broken in the vacuum and the expectation values of the meson 
fields are (cr) = /„- and (7?) = 0. With our model parameters the couplings are 
A = 20 and g = 3.3. 

Let us consider a spatially uniform system in thermodynamical equilibrium at 
temperature T and quark chemical potential fi = /is/3, where fis is the baryon 
chemical potential. In thermodynamics the most important function is the par- 
tition function, because from it every other thermodynamical quantity can be 
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determined. In general, the grand partition function reads [[75j: 



Z = Trexp - [H - fiAf) /T 



T>qVqVaT>7r exp 



fiqy°q 



(2.3) 



where J x = Jjf dr J v <i 3 x and V is the volume of the system. H and Af are the 
Hamiltonian and the number operator. In the following we adopt the mean-field 
approximation. This consists of replacing a and n in the exponent by their expec- 
tation values. In other words, we neglect both quantum and thermal fluctuations 
of the meson fields and retain only quarks and antiquarks as quantum fields. Then, 
up to an overall normalization factor: 



exp 
exp 



VU 

vu 



J VqVq exp q [i^d^ - g(cr + ij 5 r • 7?)] q + fJLqj°q 
detp { + /i7° - g(cr + i^r ■ tt) \/t) • 



The grand canonical potential is 



fi(T,/i) 



TlnZ 
V~ 



qq j 



(2.4) 



(2.5) 



where the contribution of quarks and antiquarks follows from equation ( |2.4[) and 
is 

d 3 p 



tt q q(T,/J,) = -U q 



(2vr) : 



[E + T In 



1 + e 



(ti-E)/T 



+ T In 



1 + e 



{-H-E)/T 



(2.6) 



Here f 9 = 2N c Nf = 12 is the number of internal degrees of freedom of the quarks 
with three colors N c = 3, and E = y/p 2 + M 2 is the energy of quarks and anti- 
quarks. The constituent quark/antiquark mass, M, is defined to be: 



M 



sV + Tf 2 ) 



(2.7) 



The divergent first term in equation Q2.6|) comes from the negative energy states 
of the Dirac sea. As follows from the standard renormalization procedure [73| it 
can be partly absorbed in the coupling constant A and the constant v 2 . However, 
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a logarithmic correction from the renormalization scale remains. This term is 
neglected in the following calculations. After integrating equation ( |2.6|) by parts 
the pressure, P = — fl, due to quarks and antiquarks can be written as 



P q q(T, //) = -p- [ dp^- [n q {T, fx) + n q (T, fi)] , {2.1 



COO rrA 

6^ Jo 

where n q and n q are the quark and antiquark occupation numbers, 

1 



n i\ T >y) i + e (E-n)/r 

rig(T,fi) = n q (T,-fi) . (2.9) 

The baryon-chemical potential is determined by the net baryon density 

n B = -~ = ^- 2 j p 2 dp[n q (T^) - n q (T,fi)] . (2.10) 

The net quark density is n = 3n#. The values for the a and 7? fields and thereby 
the quark masses in equation ( |2.7| ) are obtained by minimizing the potential Q 
with respect to a and jf. The results are equations of motion for the meson fields: 

dfl 

— = \( a 2 + ^- v 2 ) (r -H + gp s = 0, (2.11) 

and 

dfl 

_ = A( ( x 2 + 7f 2 -t; 2 )7f+( ? p ps = . (2.12) 

The scalar and pseudoscalar densities of quarks and antiquarks can be expressed 
0: 



as 



Ps = (qq) = gov q J -^^^[n q (T, fi) + n q (T, fi)} , (2.13) 
f d^p 1 

P P s = (qhsrq) = gnu q J 7^p^M T > P) + n g( T ^)} ■ ( 2 - 14 ) 

These densities generate the source terms in the equations of motions for the meson 
fields ( |2.11| ), ( |2.12| ). The minima of fl defined by eqs. ( 2.11|) , ( 2.12|) correspond to 
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stable or metastable states of matter in thermodynamical equilibrium where the 
pressure is maximal P ma x = —^min- 



This version of the sigma model was used earlier in [59| for thermodynamical 
calculations at nonzero T and /x = 0, and at nonzero \i and T = 0. Some useful 
formulae for the case of a small quark mass are given in Appendix [A|. 



2.3 Phase Diagram 

In our model everything is determined when the gap equation is solved 

in the (T, \x) plane. Below we present results of our numerical calculations. As 
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Figure 2.1: The phase diagram in the (//, T) plane. The middle curve is the critical line 
and the outer lines are the lower and upper spinodal lines. C is the critical point. 
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our calculations show, the model exhibits a first order chiral phase transition at 



T < T c and nonzero chemical potential. In figure |2.1| we present the resulting 
phase diagram in the (T, jj) plane. The middle line corresponds to states where 
two phases co-exist in the first order phase transition. Along this line the ther- 
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Figure 2.2: The thermodynamical potential Q, for \i = 0. The levels correspond to 
(starting from the top): T = [0, 100, 135, 155, 175, 190] MeV 



modynamical potential Q has two minima of equal depth separated by a potential 
barrier whose height grows towards lower temperatures. At the critical point C 
the barrier disappears and the transition is of second order. The other lines in fig- 
ure |2.1| are spinodal lines which constrain the regions of spinodal instability where 
(driB I 'dfi) T < 0. Information about the timescales of this instability is important 
for dynamical simulations of, for example, disoriented chiral condensate formation 



or quark-antiquark droplet formation [29 . 
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It is instructive to plot the thermodynamic potential as a function of the order 
parameter for various values of T at \i — 0, and for various values of // at T — 0. 
The first case is shown in figure TH. One clearly sees the smooth crossover of the 
symmetry breaking pattern. The energy difference between the global minimum 




Figure 2.3: The thermodynamical potentials f2 for T = 0. The levels correspond to 
(starting from the top): // = [0, 225, 279, 306, 322, 345, 375] MeV. 

and the local maximum of the potential in vacuum, otherwise known as the bag 
constant, is about 60 MeV/fm 3 . The temperature corresponding to the crossover 
transition is about 140 MeV. In figure |2.3| the same plot is shown for T = and 
a nonzero /i. Here, one clearly observes the pattern characteristic for a first order 
phase transition: two minima corresponding to phases of restored and broken 
symmetry separated by a potential barrier. It now follows that somewhere in 
between these two extremes, for some fi c and T c , there exists a critical point for 
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a second order phase transition. Indeed, this point is found and shown in figure 
|2.1| . The corresponding values are (T c , fi c ) ~ (99, 207) MeV. The behavior of the 
thermodynamic potential at fi = fi c and various T is shown in figure 2~ "A. One can 
see that the potential has only one minimum which is flattest at the critical point. 




CONSTITUENT ^ASS [MeV] 



Figure 2.4: The thermodynamical potentials f2 for fi = 207 MeV. The levels correspond 
to (starting from the top): T = [0, 50, 75, 100, 125, 150] MeV. 
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2.4 Effective Masses 



The masses of the sigma meson and pion are determined by the curvature of the 
potential at the global minimum: 



Mi 



d 2 Q 



and Ml 



d 2 Q 



da 2 ^ d7r 2 ' 

Explicitly they are given by the expressions 



1,2,3 . 



Ml 



2 a. ^ f 2 



+ 9 



2 



2tt 2 
M 2 



r 



dpp 2 



P 



+ 



1 



1 1 

+ 



and 



Mi 



TE 2 \2{l + cosh[(£ + /i)/T]) 2(1 + cosh[(£ - //)/T]) 



+ 9 



2 "9 



2tt 2 







E 



+ 



1 



1 _|_ e (E+/i)/T I + e (E-n)/T 



(2.15) 



(2.16) 



(2.17) 



Since we have set the expectation value of the pion field to zero, (7?) = 0, the quark 
mass is given by 



Mi 



2 2 

9 



(2.18) 



Figure 2J5 shows the constituent quark mass as a function of T and /1. This 
plot, of course, shows the same phase structure as discussed in the previous section. 
At jj, — the quark mass falls smoothly from its vacuum value and approaches zero 
as T goes to infinity. One could define a crossover temperature as corresponding 
to a steepest descent region in the variation of M. This again gives a temperature 
of about 140 — 150 MeV. At T = and nonzero /x the constituent quark mass 
shows a discontinuous behavior reflecting a first order chiral transition. 
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Figure 2.5: The constituent quark/ ant iquark mass as a function of ji and T. 
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The sigma and pion masses for various T and /x are shown in figures |2.6| and 
[2.7| . One can see that the sigma mass first decreases smoothly, then rebounds and 
grows again at high T. The pion mass does not change much at temperatures 
below T c but then increases rapidly, approaching the sigma mass and signaling the 
restoration of chiral symmetry. At large T the masses grow linearly with T (see 




. . . . ! . . . . L . . . . ! . . . . ! . . . . 

50 100 150 200 250 

TEMPERATURE [MeV] 

Figure 2.6: The sigma mass (solid line) and pion mass (dashed line) as functions of 
temperature for \x = (right pair) and for \i = \i c ( left pair). 

Appendix |A|) . The \i = \i c case is especially interesting. Since the sigma field is the 
order parameter of the chiral phase transition, its mass must vanish at the critical 
point for a second order phase transition, M q ~ o ~ (qq)- This means that Q has 
zero curvature at this point. Figure |2l| indeed shows 

In Fig. [2.7| the masses are plotted as function of // for T = and T = T c . For 
T = one clearly sees discontinuities in the behavior of the masses, characteristic 
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CHEMICAL POTENTIAL [MeV] 



Figure 2.7: The sigma mass (solid line) and pion mass (dashed line) in the as functions 
of chemical potential for T = (right pair) and for T = T C (left pair). 



for the first order phase transition. 

An interesting point is that the quark mass assumes its vacuum value all the way 
up to the chiral transition, and then drops to a small value in the phase where chiral 
symmetry is restored (see figure |2.5| ). This behavior is related to the appearance 
of a bound state at zero pressure. Within the linear sigma model this bound state 
was found by Lee and Wick a long time ago [^2| . Recently, it was shown in [41] 
that a similar bound state appears also in the NJL model. This behavior, however, 
depends on the value of the coupling constant g. In general, if the coupling constant 
is sufficiently large, the attractive force between the constituent quarks becomes 
large enough to counterbalance the Fermi pressure, thus giving rise to a bound 
state. To demonstrate this we have varied the coupling constant g for the sigma 
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Figure 2.8: The constituent quark/antiquark mass as a function of T for \x = (upper) 
and as a function of jjl (lower) for T = 0. The solid line represents the mass for g = 4.5, 
the dashed-dotted line for g = 3.3 and the dashed line for g = 2.8. 
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model within reasonable limits. The results are shown in figure |2.8|. It is seen in 



figure p.8| that, indeed, one can change the smooth crossover for p = into a first 
order transition by increasing the coupling constant (left panel) and change the 
first order transition in the case of T = into a smooth crossover (right panel). 
In this way a heavy quark phase comes into existence as the coupling constant is 
decreased and the bound state disappears. 



2.5 Adiabats 

In hydrodynamical simulations the entropy per baryon, S/A, is a fundamental 
quantity. One can calculate it using standard thermodynamic relations, 

S e + p-pn 

A Tn K J 

where e, p, n are the energy density, pressure and net density of quarks and 
antiquarks (n = 3ns), respectively. By studying this quantity, one can check if 
there is a tendency towards convergence of the adiabats towards the critical point 
as was claimed in |HJ|. If this was the case it would be easy to actually hit or go 



close to this point in a hydrodynamical evolution. 

Figure [2l] shows the contours of S/A in the (T, p) plane. We actually observe 
a trend which is quite opposite to this expectation. It turns out that the adiabats 
turn away from the critical point when they hit the first order transition line and 
bend towards the critical point only when they come from the smooth crossover 
region. This is understood as follows. First, note that all adiabats terminate at 
zero temperature and p = M vac , i.e. the (T,p) combination corresponding to the 
vacuum. The reason is that as T — > 0, by the third law of thermodynamics also 
S — > 0, therefore, for fixed S/A we have to require that n — > is fulfilled when 
p = M vac . For our choice of parameters the point (T,p) = (0,M mc ) is also the 
endpoint of the phase transition curve at T = 0, since the phase transition connects 
the vacuum directly with the phase of restored chiral symmetry. This is confirmed 
by figures 1213 and 



This behavior is quite opposite to the case underlying the claim in [31], where 
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Figure 2.9: The entropy per baryon number S/A. The curves correspond to (from left) 
S/A=[28,21, 17,13,11,9,6,2]. 



the hadronization of a large number of quark and gluon degrees of freedom into 
relatively few pion degrees of freedom leads to the release of latent heat and con- 
sequently to a reheating (increase of T) through the phase transition. Remember, 
however, that in our case there is actually no change in the number of degrees of 
freedom in the two phases. The change is only in the quark masses. Consequently, 
there is no "focusing" effect in the linear sigma model with quarks. On the other 
hand, the behavior of the adiabats in figure |2.5| is quite typical for an ordinary 
liquid-gas phase transition. Here liquid and gas are represented, respectively, by 
chirally symmetric and broken phases. This analogy was further elaborated in |4T] . 
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2.6 Summary 



In this chapter the results of our investigations of the thermodynamics of the chiral 
phase transition within the linear sigma model coupled to quarks were presented. 
For small bare quark masses we have found a smooth crossover for nonzero temper- 
ature and zero chemical potential and a first order transition for zero temperature 
and nonzero chemical potential. The first order phase transition line in the (T, /i) 
plane ends in the expected critical point. It has been found that the a mass is 
zero at the critical point. The behavior of the adiabats shows a pattern opposite 
to the expectations for the chiral/confinement transition of [|3T|. In fact, the phase 
transition found in this model turned out to be of the liquid-gas type. However, 
the strength of this transition depends sensitively on the coupling constants of the 
models. 

Certainly it will be interesting to use this model in hydrodynamical simula- 
tions in order to confirm or disconfirm possible observable signatures of the phase 
transition discussed in the introduction. In particular, since the model contains 
dynamical sigma and pion fields, it would be suitable to study the long wavelength 
enhancement of the fields at the critical point. 
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Chapter 3 



Evolution of Disoriented Chiral 
Condensates 

In this chapter, I review the concept of disoriented chiral condensates and the 
research done so far on this subject. Then I derive coarse-grained equations of 
motion for out-of-equilibrium condensate fields. Based on linear response theory 
the fluctuations of the thermal bath, as a response to the presence of the non- 
thermal condensate, are analyzed. Multiple aspects of the response functions are 
emphasized, setting the ground work for the chapters that follow. 

3.1 Review of Previous Results 

Formation of condensates as a result of a phase transition between a high temper- 
ature disordered (symmetric) and a low temperature ordered (symmetry broken) 
phase is common in different areas of physics. As the disordered phase is cooled, 
correlation occurs on a microscopic level, reaching a macroscopic size at the criti- 
cal temperature. In the region of the phase transition small scale physics becomes 
irrelevant and phase ordering is determined by long wavelength fluctuations. The 
time-scales involved in the transition play a major role in determining the na- 
ture of the phase ordering. When the cooling is much slower than the time taken 
for microscopic rearrangement (relaxation time), the system is said to be in local 
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thermal equilibrium, and the transition is an equilibrium phase transition. The 
QCD phase transition in the Early Universe qualifies as such, since the cooling 
rate around the critical temperature is about 1CU 5 seconds, whereas the relaxation 
time estimated from collisional processes of quarks and gluons has a very different 
order of magnitude, on the order of 10~ 23 seconds. There is also the possibility for 
the expansion and cooling to be very rapid. The long wavelengths, not adjusting 
to the fast temperature change, fall out of equilibrium. Because the time-scales 
involved in a heavy ion collision (10~ 22 — 1CU 23 seconds at RHIC, for example) 
are comparable to the relaxation time, it is now widely believed that the QGP 
confinement /chiral phase transition happens out of equilibrium. 

In a recent work [3JJ] Boyanovsky et al. discuss in a thorough and compre- 
hensive manner the relaxation of the order parameter and its phenomenological 
consequences in ultrarelativistic heavy ion collisions and in Early Universe cosmol- 
ogy. The authors focus on the possibility that long wavelength fluctuations are 
critically slowed down. With decreasing temperature the relaxation time of the 
long wavelength modes of the order parameter become extremely long. Then these 
can freeze-out even before reaching the critical point. Such a transition is said to 
be quenched. There is speculation that such phenomena can be the source of 
non-equilibrium effects in cosmological phase transitions, which may for example, 
explain the size and distribution of primordial black holes. 

The dynamics of phase ordering is described in terms of domains that form 
after the temperature has rapidly dropped below its critical value. In heavy ion 
collisions, after the plasma expands and cools through the non-equilibrium chiral 
phase transition, such domain formation is expected. These are referred to as dis- 
oriented chiral condensates, named shortly as DCC A DCC is a good candidate 
for the chiral phase transition to be directly probed. Inside a DCC domain the 
order parameter is misaligned, meaning that the chiral condensate points in a dif- 
ferent direction than in vacuum. The possible physical picture for DCC formation 
is presented in figure [D] and described in the caption. It has been suggested [[|^ 
that regions in which the chiral condensate is disoriented, corresponding to long 
wavelength fields, may spontaneously grow and subsequently decay via pion emis- 
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sion. This can be detected through the different pion decay channels: it — > 27 and 
7i^~ — > pirVp. DCC decay can lead to experimentally observable anomalies. One 
is the ratio of neutral to charged soft pions, r = — t 2 — . Each misaligned region 

o f ' no+n ch & & 

gives pions with a given r, but r is different for different regions. The probability 
distribution is P(r) = which is different than what is naively expected (|). 
Peculiar events in which such anomalies are present have been observed in cosmic 
ray experiments |46j, and are known as Centauros ||47]. Ultrarelativistic heavy 
ion collisions provide a fertile ground for studying fluctuations in r. The STAR 
detector at RHIC p5f searches for such nonstatistical or dynamical fluctuations on 
an event-by-event basis. 

The formation and decay of DCCs have been studied extensively in the last 
decade since the pioneering work of Anselm ||9| in 1989. However, it was 
first Bjorken M i arid then Blaizot and Krzywicki [||lj, who in 1992 made the 
idea widely popular. In 1993, Bjorken, Kowalski and Taylor proposed a "baked- 
Alaska" model |52] which could explain Centauro events. It was Rajagopal and 
Wilczek ||45|| , though, who first used the idea in quark-gluon-plasma physics. In 
relativistic heavy ion collisions the multiplicity of produced particles is very high. 
Due to their small mass, the majority of these are pions. Rajagopal and Wilczek 
proposed that some of these pions are the result of the coherent decay of classical 
pion fields present in these collisions. These pions would have small momenta and 
would show large fluctuations in the ratio of produced neutral pions compared 
to the charged ones. The authors assumed a perfectly quenched scenario for the 
transition: after the critical temperature is reached, the potential is assumed to 
take instantaneously its vacuum shape, and the long- wavelength fields then evolve 
according to zero temperature equations of motion. This original idea stirred much 



interest and led to the investigation of various aspects of DCC dynamics ||53||-|61 



In []62| Asakawa, Huang and Wang pointed out that quenching cannot be nat- 
urally achieved through expansion. Krzywicki and Serreau in a more recent paper 
)3| predicted a small probability for an observable DCC signal. However, an- 



other recent publication by Scavenius and Dumitru [|30| reports the realization of 
quenched initial conditions when the chiral phase transition is of first order. 
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Figure 3.1: TOP: The T = vacuum potential has the shape of the bottom of a wine 
bottle and the ground state is degenerate; MIDDLE: For T > the hill becomes more 
and more shallow; BOTTOM: At T > T c symmetry is restored and there is a unique 
ground state. It is assumed that in a heavy ion collision such a state has been achieved. 
Due to rapid expansion T < T c potential is TOP again. The symmetric state becomes 
unstable and will roll down to a stable ground state. This roll down in radial direction 
can develop configurations in azimuthal direction. These are DCCs. 
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A quench in the strict sense defined by Rajagopal and Wilczek, is a very dras- 
tic assumption which leads to a complete decoupling of soft and hard modes. In 
a relativistic heavy ion collision, though, pions are produced in a very wide mo- 
mentum range. Soft pions, with small momentum, are numerous, but carry little 
energy. Hard modes, having momenta of a few hundreds of MeV, carry most of the 
energy. Their effect should certainly not be neglected. Biro and Greiner |56| stud- 
ied first the effect of hard modes on the non-equilibrium evolution of soft modes. 
They proposed hard pions as a thermalised background and considered the chiral 
condensate in this background. The presence of a heat bath leads to dissipation 
that shows up in the equations for the condensate fields. The result of is that 
dissipation suppresses the DCC formation. The problem with their model, though, 
is that they calculated the dissipation coefficient in </> 4 theory based on p7|, but 
this is good only in the symmetric phase. In the broken phase the masses are 
not degenerate anymore and also the 3 interaction is allowed. In | 68fl , a work by 
Rischke, these problems have been eliminated. He accounted for cubic interaction 
and so for the possibility of decay, but neglected collisional interactions, arguing 



that these are higher order in the coupling. Csernai, Ellis, Jeon and Kapusta |70 
showed that scattering processes can be significant for dissipation. 

All of the above mentioned studies looked at homogeneous condensates using 
different approaches: in Langevin-type equations of motion; in [^HJ closed time 
path method; in |)8| influence functional formalism; in [[7(J linear response theory. 
The relaxational dynamics of inhomogeneous condensates have been investigated at 
one-loop level by Boyanovsky, D'Attanasio, de Vega and Holman The authors 
concentrate on the evolution of the sigma component only, but they assume the 
vacuum value for the order parameter. We know that the equilibrium value of 
the condensate is changing with temperature, and their assumption is correct only 
at low temperatures. Moreover, in their model the mass of the pion is given by 
hand, so this is not a consequence of symmetry breaking. The effect of the heat 
bath on the low-frequency modes is discussed at one-loop level by Jakovac, Patkos, 
Petreczky and Szep p9| too. 



The work presented in this dissertation was substantially influenced by |70 
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where the dynamical evolution of the sigma condensate has been followed within 



the linear sigma model. My analysis is an extension of [70] allowing for DCC for- 
mation. I present a microscopic, field theoretical description of the non-equilibrium 
dynamics of soft homogeneous and inhomogeneous sigma and pion field configu- 
rations coupled to a thermal bath. To the best of my knowledge, a consistent 
incorporation of relaxation processes at one- and two-loop level has not been done 
before in this context. In this model dissipation follows from first principles. This 
means that the dissipative terms in the equation of motion arise directly from the 
Lagrangian. Moreover, the temperature-dependence of the order parameter natu- 
rally follows, too. A convenient framework for my DCC studies is the linear sigma 
model. It is convenient because: 

• it is formulated in terms of the chiral field = (a, 7?); 

• it exhibits symmetry breaking such that in the ground state the scalar con- 
densate is nonzero, o ~ (qq) 7^ 0, and the pseudoscalar condensate is zero, 
tT ~ (qf^q) = 0; 

• it implements a phase transition at finite T c such that for T > T c chiral 
symmetric is restored and (qq) = 0; 

• it is a renormalizable theory. 



3.2 Equations of Motion 

The dynamics of the condensate are completely determined by the evolution equa- 
tions in space and time for the long-wavelength fields. In the following, the coarse- 
grained equations of motions are derived. First, the theory with exact chiral sym- 
metry is analyzed. This is an idealized case, which significantly simplifies some of 
our formulas, and allows for the study of massless bosonic waves, the Goldstone 
modes. Such modes are present not only in particle physics, but also in condensed 
matter physics (e.g. spin waves). The low temperature behavior of the pion is well- 



known from chiral perturbation theory [12| and, accordingly, the pion is massive 
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and chiral symmetry is explicitly broken. Therefore, for a completely consistent 
discussion one should work with massive pions. This is the topic of the second 
section. 



3.2.1 Chiral Limit 

The Lagrangian of the O(N) sigma model with no explicit symmetry breaking term 
is given by 

L(a,n) = 1 -(da) 2 + ±(dny-U(a,n) , (3.1) 

where a and 7r are the sigma and the N — l component pion fields, and the potential 
is given by 

U(a^) = \^ + ^-flf . (3.2) 



Note that the fermionic part of the original sigma model Lagrangian ( |1.12|) has 
been neglected here. This is allowed since our focus is on the condensate formed 
in the symmetry broken phase where quark degrees of freedom are confined al- 
ready. In the hadronic phase a natural choice would be nucleons. However, we 
justify their exclusion by their mass, which is much larger than the other scales 
involved. Therefore, processes that involve nucleons would be strongly suppressed 
with respect to pure bosonic interactions. The evolution of the fields under ther- 
mal equilibrium circumstances has been well studied by several authors (see for 
example WW). The usual procedure is to express the fields as 



a(t, x) — v + (Tf(t, x) and 7T;(£, x) = iTif(t, x), i — l,...,N — l . (3.3) 

Here v is the chiral order parameter, the thermal average at temperature T, chosen 
to lie along the sigma direction. Of is the fluctuation about v, along the radial 
direction and Hif is the fluctuation of the ith component of the pion field about 
zero, in the azimuthal direction. 

As explained in section |3.1] , there is good reason to believe that a more realistic 
picture involves an out-of-equilibrium evolution of the fields. The approach we 
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take is a semiclassical description. The idea is that instead of (|3.3|), we decouple 
the fields into their low and high frequency modes. This can be done, for example, 
by introducing a momentum scale A c , as in [J7IJ. Then one integrates out the high 
frequency modes, and obtains an effective classical theory for the long wavelengths. 
The fields can be separated as follows: 

a(x) = a(x) + (Tf(x) 

TTi(x) = 7Ti(x) + 7Tif(x) , Z=l,...,iV-l (3.4) 

• a and tt are slowly varying condensate fields, representing low frequency 
modes with momentum | k |< A c , occupied by a large number of particles. 
These soft modes lie entirely in the non-perturbative regime, and may be 
treated as classical fields. 

• Of and mf are high frequency, fast modes, with | k |> A c . These hard modes, 
representing quantum and thermal fluctuations, constitute a heat bath. 

The problem to be solved, then, is to describe the evolution of low frequency 
classical fields that are embedded into a thermal bath. In our approach, the soft 
modes follow classical equations of motion, whereas the effect of the hard thermal 
modes is taken into account in a perturbative manner. We are well aware that 
any perturbative treatment can be a target of criticism, since the sigma model is 
a strongly coupled effective theory. However, one can still determine qualitative 
aspects using this model. 

Euler-Lagrange field equations are derived after inserting ( |3.4j ) in the Lagrangian 
( |3.1| ). For the sigma field one obtains 

N-l N-l N-l 

d 2 a + A5- 3 + 3\ap + \a tt| + S\a f a 2 + \a £ Trfy + 2\a £ ^hif ~ ^& 

i=l i=l i=l 

N-l N-l N-l 

+ AaJ + \a f ^ + Act/ ^ nf f + 2\a f ^ nfr if - \f 2 w a f = , (3.5) 

i=l i=l i=l 

and for the ith component of the pion field 

N-l N-l N-l 

j^i=l j^*=i i^»=i 
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N-l N-l N-l 

+ Xltff + X-Kif Kj+^if Y 7r| / + 2A7T i / Y TTj^/ + A7TiCT 2 + Ao- 2 ^/ 

j^i=i j^»=i i^ i= i 

+ \ajfci + Xajn if + 2Xaafn i + 2Xaaf'n- i f - Xf^i - Xf^TTif = . (3.6) 

Now we average the equations of motion over time and length scales that are short 
compared to the scales characterizing the change in the slow fields, but long relative 
to the scales of the quantum fluctuations. This is known as coarse-graining. The 
thermal average of high frequency fluctuations is (07(2;)) = and (irf(x)) = 0, 
while (cr(x)) = a(x) and (7ti(x)) = 7Tj(x). Throughout this work the notation (O) 
refers to the non-equilibrium ensemble average of an operator O. The equilibrium 
ensemble average is denoted by (0) eq . It should be noted that we allow for a 
nonzero ensemble average not only along the sigma direction, but also in the pion 
direction. In other words, we allow the formation of disoriented chiral condensates. 
The coarse-grained equations of motion are 

N-l N-l 

d 2 a + Xa 3 + 3X(a 2 )a + Xa Y tt 2 + Xa Y (n 2 f ) - Xf 2 a 

i=i i=i 

N-l 

+2A]T( ( 7 / 7r, / >^-A/ 2 a = , (3.7) 
i=i 

and 

N-l N-l N-l 

9 2 ^ + A^ 3 + 3A(7r 2 / )^ + A^ Y ^ 2 + A^ E (4> + 2A £ nA^jf) 

j^i=l i^*=i j¥»=i 

+A^a 2 + A7f i <7r? / > + 2Xa(n if a f ) - Xj% = . (3.8) 

The cross correlations between the fluctuations of different fields are nonzero, 
(<Pi<Pj) 7^ 0; however, when averaging, all cubic fluctuations of the form (ipifjipk) 
were considered zero, meaning that only one-loop order is included. Later in this 
section the extension to two-loops is presented. 

In the broken symmetry phase, T < T c , we can separate the non-equilibrium 
condensate fields into their equilibrium value and a slow fluctuation about this. 
Choose the equilibrium condensate along the sigma direction: 

a(x) = v(x) + a s (x) 

TTi(x) = ir is (x), i = 1, ...,N - 1 . (3.9) 
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For now, let us assume that this equilibrium condensate is static and homogeneous, 
v(x) = v. Clearly, this need not be the case, and we will address this issue 
later. The thermal equilibrium ensemble of the hard fluctuations is affected by 
the presence of the condensate. Full ensemble averages, which are basically the 
two-point functions of the thermalized fields, can be written as 

<*?> = + 8(a)) 

faff) = faff) eq + S faff) 
(a f n if ) = 5(a f n if ) 

(KifKjf) = Sfaif^jf) ■ (3-10) 

The deviations in the fluctuations, 8(a"j), 6 faff), Sfaf^f), and SfaifiTjf), are the 
responses of the fast modes to the presence of slow a s and n s background fields. 
These responses are proportional to the slow modes raised to some positive power. 
A detailed discussion, including an explicit evaluation using linear response the- 



ory, is presented in section |373| . Notice the absence of faf^if) e q and faifXjf)eq- 
The reason is that the correlation functions in equilibrium, in the absence of a 
background, are diagonal. 

Equations (|3.7| ) and (|3.8|) must also hold in equilibrium, at fixed temperature. 
Setting a s = and n s = results in an equation satisfied by the equilibrium 
condensate: 

JV-l 

\v 3 + 3\v(a 2 f ) eq + \vJ2 (4)e 9 - A/> = . (3.11) 

i=i 

This has two solutions: 

v = , T>T C , (3.12) 

and 

,2 p2 q / \ f AT 1 \ 



V 



fl - 3(aj) eq -(N- l)fai f ) eq , T<T C . (3.13) 



We are interested in the second solution ( |3.13|) , since this represents the low tem- 



perature, symmetry broken phase. Since all components of the pion field are equiv- 
alent, we could write the sum as Sill 1 faff) eg = fa})eq = {N — l)faff) eq , where 
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now TTij is one of the components. This is a convention that we adopt throughout 
the rest of the paper. The equations describing the evolution of the condensate 
fields are readily obtained after inserting ( |3.1(J| ) and ( |3.9| ) into ( |3.7D and 



d 2 a s + 2Xv 2 a s + Xv 35 (erf) + (N — 1)6 (nf) = -X (a 6 s + 3va 2 s + (N - l)vn 



+ {N - l)va s 7r 2 + 2(N - l)S(a f 7r f )7r s + 3a s 5(a 2 f ) + (N - l)a s 5(vrj)) (3.14) 



and 



d 2 rc is + 2Xn is (n 2 f ) eq - (<Tf)eq + 2\v5{iTif(Tf) = -X (vrf s 
+ {N - 2)7c is 7r 2 + 2a s 8(ir if <j f ) + 3ir is 5(Trf f ) + 7c is 5(a 2 f ) 
+(N - 2)n is 5(7i 2 f ) + 2(N - 2)7r s 6(n lf n f ) 



2va s Tr is + a s 7T is 



(3.15) 

The above set of equations is a coupled non-linear integro-differential equations (see 
section |3~3|) , and it is hardly possible to solve, even with access to supercomputers. 
Let us assume only small deviation from equilibrium, |<r s |, |7r s | v, neglecting 
higher order terms. The linearized field equations read 



d 2 a s + 2Xv 2 a s + Xv 35(a 2 f ) + (N - I)S(ttJ) 







and 



d 2 n s + 2Xn s (n 2 f ) eq - (aj) eq + 2Xv5(n f a f ) = 



(3.16) 



(3.17) 



Remember, ir s and itf refer to one of the N — 1 components of the pion field. The 
three physical pions, ir° and 7r =t , are the degrees of freedom in the N = 4 model. 
Also, the above results were obtained with the assumption of fixed temperature, 
and so the equilibrium condensate was constant. In reality, the temperature drops 
as time elapses, and v changes with it (see previous discussion in section |3~T| ). 

We close this section by stating the results when accuracy is extended up to 
two-loop order. Technically this means that, when averaging equations ( |3.5| ) and 
(|3.6j ), we need to keep the three-point functions too. The ensemble average of the 
product of three fields is given by: 



(a}) = (a)) eq + 8(a)) 



38 



(a f n 2 f ) = (a f ir 2 f ) eq + 5(a f n}} 
W> = 8(1$) 

(a 2 f ir f ) = 5(a 2 f 7r f ) . (3.18) 

On account of cubic couplings one finds nonzero (crj) eq and (afir 2 ) eq . These are 
evaluated in the usual way (see section |3.3| ) by performing an expansion in the 
coupling constant of the exponential of the action |73| . When including these 
terms, the equilibrium value of the condensate is modified compared to its previous 
expression (|3.13|): 



v 2 = fl ~ S(aj) eq — (N — l)(nj) eq - _ (iV _ x) <W>3 . (319) 

J v v 

The linearized equations of motion acquire some extra terms compared to equations 

( CT ) and Q3T7D : 



d' 2 a s + rfi 2 a a s + \v 35 (a j) + (N - l)6(rf) + \S{aj) + \{N - l)5(a f n 2 f ) = 



(3.20) 



where 



and 



where 



ml = x[2v 2 - { -^-(N-l) { -^I^) , (3.21) 



d 2 7T s + 2Xv5(a f 7T f ) + m> s + \5(tt}} + \S(ap f ) = , (3.22) 



K = X (2tfU - 2(a 2 f ) eq - { ^-(N- . (3.23) 

\ 1 v v ) 

3.2.2 Explicit Chiral Symmetry Breaking 

The explicit breaking of symmetry is implemented in the Lagrangian through the 
presence of a term linear in the sigma field (see discussion in |1.2| ). The potential 
term is 

U(a, n) = j (a 2 + t? 2 - v 2 ) + Ha . (3.24) 
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Since the procedure to obtain the effective equations of motions is identical to the 
above derived massless pion case, we are going to sketch the major steps of the 
derivation, giving details only where the physics is different: i) separate the low 
and high frequency modes of the fields as in ( |3.4|) ; ii) write down Euler-Lagrange 
equations; iii) average these with respect to the statistics of the hard modes; iv) 
separate the condensate fields into their equilibrium value and the slow fluctuations 
about these as in ( |3.9|) ; v) break the hard fluctuations into a sum of equilibrium 
fluctuations and deviations from these as in ( p.lO|) ; vi) drop all terms other than 
linear in the fields. 

One of the significant differences between the two cases is in the behavior of 
the equilibrium condensate, which now at one-loop order satisfies the following 
equation: 



Xv 3 + X (3(a 2 f ) eq + (N- l)(ir}) eq ) v - Xv 2 v 



H = 



(3.25) 



Similarly, notice, that unlike before, v — is not a solution of ( |3.25| ). This means 
that the condensate is never fully dissolved by fluctuations, not even at high tem- 
peratures, and the symmetry is never restored. This should come as no surprise, 
though, since including a finite pion mass means that the O(N) symmetry is only 
approximate. The linearized equations of motion are 



d 2 a s + mla a 



Xv [3S(a 2 f ) + (N - 1)5(ttJ) 



= 



(3.26) 



2Xv 2 + — 



(3.27) 



and 



d 2 ii s + m 2 TT s + 2Xv5{a fit f) + 2A7r s {it 2 ) e q — 



f/eq 







(3.28) 



- 2 H fn 2 

m„ = — = — 



V V 

When the strength of the symmetry breaking diminishes, H 
the previous section are recovered. 



(3.29) 
0, the results of 
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3.2.3 Discussion 



Equations ( |3.16| ) and ( |3.17| ), derived above for the theory with exact chiral sym- 
metry, and ( |3.26| ) and ( |3.28|) when the symmetry is explicitly broken, constitute 



the basis of the study presented in the rest of this dissertation. Let us enumerate 
some of the focal points: 

• The mass of the sigma meson, m a , is given in terms of the equilibrium con- 
densate, v , and as such, it is temperature dependent, since v includes thermal 
fluctuations of the fields. However, the fluctuations themselves are depen- 
dent on the sigma mass. Therefore, in order to determine the temperature 
dependence of the mass, one has to perform a self-consistent evaluation. This 
analysis is presented in Chapter f|. 

• In the theory with exact chiral symmetry the mass of the pion is zero. Gold- 
stone's Theorem requires that it remains zero at every order in perturbation 
theory. However, equation (|3.17| ) for the pion condensate shows a mass, 



mi = 2A 



(7rj) e5 — {(T 2 f) eq 7^ 0, that is non- vanishing at finite temperature. 
This is an apparent violation of Goldstone's Theorem. It is apparent only, 
because in Chapter [| it is shown that Goldstone's Theorem is not violated. 

The effect of hard modes on the condensate enters through the deviation of 
the field fluctuations from their equilibrium values. This is nothing but the 
response of the fast modes to the presence of the slow ones. In this work 
these deviations are identified as time-delayed response functions, and are 
evaluated using linear response theory. The response functions renormalize 
the equations of motion and give rise to dissipation. The consequences of 
such renormalization is the modification of the properties of particles in the 
condensate, such as mass, and changes in their speed of travel (see Chapter 
Due to the possible interaction with the hard modes decay channels 
open up and particles can scatter. These processes are responsible for the 
dissipation of the condensate. For detailed discussion see Chapters |j] and 0. 
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3.3 Response Functions 



Studying the dynamics of the condensates means solving equations Q3.16Q and 
(13.171 ) in the chiral limit, and equations (|3.20|) and (|3.22| ) for the theory with 
massive pions, together with some initial conditions. An important step towards 
this goal is the evaluation of the response functions. 

Linear response theory is a very convenient tool when one is interested in 
monitoring the effect of an an external field applied to a system initially in thermal 
equilibrium and this effect is small (for a comprehensive introduction see |[74|| ). 
The response function expresses the difference between the expectation value of an 
operator before and after an external perturbation has been turned on. Consider 
a perturbation described by H per turbed, and turned on at time to- 

5{6{t)) = (6{t)) perturbed -{d(t)) 

= i dt' { H pert {t') } 0{t) ) unperturbed ■ (3.30) 

The expectation value is evaluated in the unperturbed ensemble. Extension of this 
expression to spatial coordinates is straightforward. 

Our program is to use the result of linear response theory to evaluate the effects 
of the soft condensate fields on the thermal medium. According to (|3.30|) , these 
are just the commutators of the field operators at different powers evaluated at 
two separate space-time points in the fully interacting, but unperturbed ensemble. 
The initial conditions are set by the assumption that in heavy-ion collisions the 
system reaches a state of approximate local thermal equilibrium, then it cools 
while expanding and reaches the critical temperature. We define the initial time 
t = by when the critical temperature is reached. This implies a s (0,x) = and 
7? s (0,x) = 0. The general expressions to be evaluated are: 

5(a n f {x)^{x)) =ifdH J d 3 x'([H perturbed (x'), a n f (x)nf (x)]) eq (3.31) 

The respective powers n and m are identified from (|3.16| ), ( |3.17| ), ( |3.26|) , and ( p. 28 



Let us emphasize that the responses should be evaluated in the unperturbed, equi- 
librium ensemble, which does include all the interactions between the different 
modes. 
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3.3.1 Couplings Between Different Modes 

The possible couplings between low and high frequency modes are determined by 
evaluating the potential U(a, n), with the fields separated into their slow and fast 
components: 

a = v + cr s + af , (3.32) 

7? = 7T S + 7?/ . 

The resulting Hamiltonian contains positive powers of the a s and tt s fields. For 
small departures from equilibrium it is enough to keep the dominant, at most linear 
terms only. With this assumption, the relevant couplings are the following^: 

H Usaf = X(3v 2 - fl)o s o f + 3Xva s a 2 + Xa s aj 
H^ f = X(v 2 - fl)rr s -K f + Xn s ii 3 f 
H as7rf = Xva s 7i 2 f 

H-KsOfTTf = XTx s a 2 -Kf + 2Xvn s o-fiif 

H a3 a f TT f = X(T s a f 7i 2 f (3.33) 

As before, tt Sj / refer to one of the — 1 components of the pion field. Next, 
we focus on the detailed evaluation of the two-point functions of the thermalized 
fields. According to the our treatment, these have only linear dependence on the 
non-thermal fields, and are obtained by inserting ( 3.33| ) in ( 3.31|) . The technical 



details of this derivation are uninteresting and so we only state the results: 



S(aj(x)) = 3iXv f dt' I d 3 x'a s {x')(\a 2 f (x'),a 2 f {x)\) eq . , (3.34) 



o 



5{iif(x)) = iXv / dt' d x'a s (x')( 7if(x'),7rf(x) ) 



'eq. ■> 



(3.35) 



5{a f(x)ix f(x)) = 2iXv J dt' J d x'tt s (x' ) ([a f(x')7Tf (x'), a f(x)nf(x)]) eq , . (3.36) 



1 These couplings refer to the massless pion case. The couplings involving massive pions are 
the same with the replacement f% — > f% — rn^/X. 
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The three-point functions can readily be read off when inserting ( |3.33| ) into equa- 
tion Q3.31 ). In what follows, the evaluation and physical interpretation of the 
expectation values is pursued. 



3.3.2 The Expectation Values 



Obtaining explicit expressions for ( |3.34j ), fl3.35| ), and fl3.36| ), requires the evaluation 
of the expectations values of the commutators. Based on the commutator relation 



[AB, CD] = A[B, C]D + [A, C]BD + CA[B, D] + C[A, D]B 
the results are 



(3.37) 



([a 2 (x'),a 2 (x)\) eq . 
([rf(x'),n 2 f (x)\) eq . 
([a f (x')Tc f (x'), a f (x)7r f (x)]) eq . 



2(D<(x,x') 2 -D>(x,x') 
2(D<(x,x') 2 -D>(x,x') 
D<(x,x')D<(x,x') - D>(x,x')D>(x,x') . 



Keeping in mind that t' < t, where t is the time elapsed after switching on the 
perturbation, and t' is the time-variable that has its values in the [0, t] interval (see 
( |3.30| )), the following notation has been introduced: 



D>(x,x') = (a f (x)a f (x')) eq . , D>(x,x') ee (n f (x)n f (x')) 



eq. i 



and 



D<(x,x') ee (a f (x')a f (x)) eq . , D<(x,x') = (TT f (x')n f (x)) eq . . (3.38) 



The functions D > and D < are known to define the spectral function, which is the 
inverse Fourier-transform of the thermal average of the commutator: 



p(k) 



D > (k)-D < (k) 

d A xe ik{x - x ' ] (D>(x, x') - D < (x, x')) 

dV^Mx),^x')}) eq . , 



(3.39) 



where (f = 07, iTf, and x = (t, x) and k = (k°, k) are four-vectors in coordinate and 
momentum space. This spectral density is a crucial ingredient of finite temperature 
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field theories, since it determines the real time propagator, otherwise known as the 
thermal Green's function: 

d 4 k 



D(x,x') = J ±JL e -*WD{k) 



d 4 k 



Other useful relations are 



D>(k) = (l + f(k°)) Pi (k) , 

Df(k) = f{k°)pi(k) , z = a,n , (3.41) 

where f(k°) = [e k °^ — is the Bose-Einstein distribution and (3 = T~ l . Signifi- 
cant physics is carried in ( |3.40| ): the poles of this real-time propagator correspond 
to quasi-particle (collective mode) excitations. The energies and lifetimes of the 
mesons of our interest are determined in Chapters |] and respectively. More on 
the properties of spectral densities is found in [J76| . Another approach to finite 
temperature field theory, the imaginary time approach, is presented in . 

The final results of this section are the two-point functions of thermalized fields 
in terms of spectral functions, obtained by the insertion of ( |3.41| ) and (|3.38| ) in the 
expressions ( |3.34| ), ( |3.35| ), and ( |3.36| ): 

5(a}(x)) = i3\v j'dt' J d 3 x'a s (x')2(D<(x,x') 2 - D>(x,x'f) 

= -z6Xv [ d*x'a s {x') f f J^L P -M(*-*') 



(2tt) 4 7 (2tt) 4 

xp CT (p)p CT (g)(l + /(p ) + /(g )) (3.42) 

xp.(pK(g)(l + /(p ) + /(<?°)) (3.43) 
5(a f (x)7r f (x)) = -i2\v f d A x'<K s {x') I -p- I ^- e -^)(—') 



(2tt) 4 7 (2tt) 

xp.(p)p.(g)(l + /(p°) + /(g )) • (3.44) 
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Chapter 4 
Meson Masses 



As particles propagate through matter, they become dressed by the interactions 
with the constituents of the medium. Dressed means that they acquire an effective 
mass, which is the physical mass, and which is different from the one measured 
in vacuum. In this chapter a self-consistent evaluation of the meson masses at 
finite temperature is presented in the model with and without explicit symmetry 
breaking, as introduced in chapter |3.2j . Special attention is given to the fulfillment 
of Goldstone's Theorem in the theory with exact chiral symmetry. 

4.1 Sigma Mass 



Recall the equation of motion for the classical sigma field, derived in section |3.2.1 
d 2 a s + mla s + Xv [35(aj) + (N — l)6(n 2 f )] = , (4.1) 

where 

m 2 a = 2Xv 2 , (4.2) 
is the mass of the sigma meson, given in terms of the equilibrium condensate, 



2 - fl-M)e q -{N-l){n)) eq . (4.3) 
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The sigma mass is temperature dependent since it includes the thermal fluctuations 

r 2 ) eq and (ttJ) eqi 



of the fields. These equilibrium field fluctuations, {(T 2 r) eq and (7r?) e9 , are determined 



from the general form for a two-point function 

(<p{x)<piy)U = (v(xMy)e Sl )o = J -^e- ip[x - y) D > {p) , (4.4) 



where the time-components satisfy the condition x° > y°. Here Q eq means fully 
interacting equilibrium ensemble average and ()o refers to the non- interacting, free 



ensemble average. The usual way for finding an explicit expression for (fOj) is to 
expand the interaction-action, Si, in powers of the coupling [7JJ. As a first approx- 
imation, one can keep only the leading order term of the perturbative expansion 
in A, meaning that the propagator uses the bare mass. The scalar meson field 
fluctuation is then 

<^=/(S^ (i+2/ - (B)) - <«) 

Here f a is the Bose-Einstein distribution function of sigma mesons with energy 



E = yp 2 + m 2 , and m a is the zero temperature bare mass. The first term in (|4.5|) 
is the vacuum contribution, while the second term is due to finite temperature 
effects. The zero temperature part is divergent in the ultraviolet limit. This 
divergence can and should be removed using vacuum renormalization techniques 
J73[. The finite temperature does not introduce any extra divergence since it is 
regularized by the distribution function. Therefore, T = renormalization is 



enough to obtain finite results (see also |75], page 43]). Since our focus is on 
the physics at finite temperatures, in what follows we are going to neglect zero 
temperature contributions, assuming that these have been renormalized away. 
The pion field fluctuations have a similar form: 



It is important to mention that the momentum integration in (|4.5|) and (|4.6|) has 



a lower limit, A c , due to the restriction to hard momenta, \p\ > A c . Therefore, the 
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equilibrium distribution functions have the general form of 

1 



f(E) = -g 0(p - A c ) , E = yjp* + m 2 . (4.7) 

er — 1 

So far m is the zero temperature renormalized mass. A complete treatment, how- 
ever, involves the full propagator that uses the dressed mass. Accordingly, the 
fluctuations of the thermalized fields are themselves dependent on the sigma mass. 
The temperature dependent mass, rh a , results as the self-consistent solution of ( |4.3| ) 
with (|4.5| ) and ( |4.6|) , known in the literature as the gap equation. Self-consistent 
solutions are approximate only, because in general they resum only certain classes 
of diagrams in the perturbation series. For example, self-consistent resummation 
of tadpole diagrams is known as the Hartree approximation. At this point, it 
is worth mentioning that using a self-consistent approximation makes the usual 



renormalization procedure difficult, as discussed recently in [[f7]]. The argument, 
according to which the renormalized divergent term can be neglected, because it is 
independent of temperature, is really correct only when the mass is the bare mass. 
A self-consistent calculation involves the temperature-dependent mass, leading to 
the temperature-dependence of the divergent term. Renormalization thus results 
in temperature-dependent renormalization constants, and these should not be ig- 
nored. However, such a treatment is beyond the scope of this dissertation, and in 
what follows, we are going to ignore the divergent term, while still being aware of 
this approximation. 



4.2 Pion Mass 

"If there is a continuous symmetry transformation under which the Lagrangian is 
invariant, then either the vacuum state is also invariant under the transformation, 
or there must exist spinless particles of zero mass. " This statement, written by 
Goldstone, Salam and Weinberg in |l4j], today is most commonly known as Gold- 



stone's Theorem, and has to be satisfied at every order in perturbation theory. 

The spinless, massless excitations, also called Goldstone bosons, in QCD are the 
pions, and are associated with the phenomenon of spontaneous chiral symmetry 
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breaking. In the exact chiral limit we know that pions are massless at the tree 
level, m n = 0. In the following we verify if this stays the same to first order in A. 
Let us look at the linearized field equation derived for soft pions in section |3.2.1| : 



d 2 ii s + m\^ s + 2\v5(7if<jf) 
where the notation introduced is 







(4i 



2A 
2A 



(^f)eq ~ (Of), 



d 3 p 



(2vr)3 2E a 



(1 + Va) 



(4.9) 



To fulfill the requirements of the theorem, expression ( |4.9|) should vanish at all 
temperatures. However, there is no reason for this to happen, except at zero tem- 
perature, where the thermal fluctuations themselves vanish, or when the masses of 
the pion and the sigma are equal, which is expected only at the critical tempera- 
ture. Thus, at first inspection one finds a violation of a basic theorem. This would 
mean that the presence of the heat bath alters the Goldstone boson nature of the 
pion, leading to the absence of massless excitations in the theory. This problem 
needs to be dealt with, because Goldstone's Theorem should not be violated under 
any circumstances. The same problem has been found also in J7^, BIJ and has 
been fixed using the random phase approximation in |78| |, and direct but tedious 
evaluation in |59| . In the following, we present a simple and clear way to prove 
that this violation is only apparent. 

First, one has to recognize that (|4.9| ) includes the one-loop tadpole contributions 
only. There is another one- loop diagram that contributes to order A. This is 
incorporated in the equation of motion through the response function 8(iTf<jf) (see 
section ^|), given by fl3.44| ). Since the response is already of the order of Xv ~ A 1 / 2 , 
and there is an overall factor of 2At> ~ A 1 / 2 in front of it, the expectation values 
can be evaluated at lowest order, meaning that the interacting ensemble average 
(...)eq. can be replaced by the free ensemble average (...)o- This is equivalent to the 
insertion of the free spectral function, 



Pfree(p) = 2%e(p°)5(p 2 - VT? 



(4.10) 
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into expression ( |3.44| ). After evaluating the frequency integrals: 



5(<Jf7Tf) = i2\m (J \ d x'7r s (x') 



d 3 p f d 3 q 



i(p+q)(x-x ') 



(2tt) 3 J (2tt) 3 



F(p,q,t') , (4.11) 



with 



F(p,q,lf) 



(! + /,+ U) (e t{E ° 



+E 7r )(t-t') _ -i(E„+E„)(t-i?) 



-E^)(t~t') _ e -i(E<r-E n )(t-t') 



(4.12) 



Since the deviation from equilibrium is assumed to be small, one can approximate 
the condensate with the first term from the Taylor expansion about its equilibrium 
value, tt s (x') ~ tt s (x), and can be taken out of the integral. Now, the field equation 
has a simplied form: 

d 2 7t s (x) +mlir s (x) = . 



The total pion mass is then 



where 



a\ = i2Xm 2 r I dt' 
= i2\ml / dt 
= i2\ml 



d 3 p f d 3 q 



(4.13) 



(4.14) 



o J (2n) 3 J (2tt) 3 
d3p ft dt'F(p,t') 



(2tt) 3 Jo 



(4.15) 



and Ei = yp 2 + vn? a and E% — \p |. Since the coherence time is finite, one can 
extend the time integral from zero to negative infinity^. 

rt rt rQO 

/ dt'e ±iE ^ = / dt'e ±iE ^ = / dTe ±iET 

JO J~oo Jo 

= H dTe^ E±ie)T = (4.16) 

Jo TiE + e K 1 

1 Another argument: for negative times, t' < 0, there is no condensate, tt(x') = 0, because this 
time period corresponds to the symmetric phase, above the critical temperature. 
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where e — > + is a small positive real number. Then 



a i 



i2\rh 2 



d 3 p 



2X 



°J (2ir) 3 AE a E w 

d 3 p r i 



2^(1 + 2/^-2^(1 + 2/,) 



(27T) 



2^ 



l + 2/ CT 



2-Ett 



;i + 2/.) 



Insert ( |4.17 ) and Q4.9| ) in the expression for the true pion mass, ( |4.14 ): 



2A 

+2A 




d 3 p 1 
(2tt) 3 2E n 
d 3 p 
]2nf 



• 1 + 2/ "»-/(Si2i (1 + 2/ " 



2£„ 



;i + 2/ CT 



(2tt) 3 2£ ct 
1 



2£L 



1 + 2U) 



(4.17) 



(4.18) 



thus proving that pions stay Goldstone bosons at one-loop level. Another simple 
and straightforward proof is presented in section |5.1.2| when working in frequency- 
momentum space. 



4.3 Numerical Results 
4.3.1 Chiral Limit 

In the theory with exact chiral symmetry, m n = 0, the sigma meson mass is 

ml = 2\[f*-3(a 2 f } eq -(N-l)(7r}} eq \ . (4.19) 
The hard sigma field fluctuation, 

<°?> = 5?j[ *¥3^T- (4 - 20) 



should be calculated with the meson field mass m CT , that is E = yp 2 + m 2 . For 
a massless pion field the fluctuation can be evaluated analytically, resulting in a 
simple form: 

<* ?> = • (4.21) 
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Figure 4.1: Temperature dependence of the condensate in the 0(4) model with 
exact chiral symmetry. The value of the coupling can modify the order of the 
phase transition. 
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When evaluating the field fluctuations of hard modes, the approximation A c — > 
has been made. Small momentum modes are very few inside the hard momentum 
bath, and including their contribution to the integrals introduces negligible small 
error. Physically, A c = means that the condensate fields contain only k = 
modes, meaning that the classical field configurations are homogeneous. In general, 
this needs not to be the case. For more discussion and analysis of k ^ modes see 
Chapter || on pion propagation. 

At T = thermal fluctuations obviously vanish, and the familiar w? a = 2\f% 
zero temperature mass is recovered. The coupling constant is chosen to be A = 18, 
corresponding to a vacuum sigma mass of m a = 558 MeV. The order parameter, 
which is the equilibrium condensate v, not only determines the mass of the sigma 
but also the order of the phase transition. The behavior of the condensate as the 
temperature changes, determined numerically in a self-consistent manner, is rep- 
resented by the solid line in figure [4.1| . In this plot N = 4. The effect of thermal 
fluctuations is to decrease the condensate with temperature from its vacuum value 
of 93 MeV. At a critical temperature of T c ~ 130 MeV the condensate is completely 
dissolved. Above T c the fluctuations are much too large to allow the formation of 
any condensate. The critical temperature, defined by the vanishing of the conden- 
sate, can be estimated from the expression (|4.22| ), and is given by T c 2 = j^r^f^- 



In the 0(4) model T c = V^f^ = 130 MeV, which coincides with our numerically 
determined value. Note also the positive-definite mass at all temperatures. The 



tachyon problem present in mean field approximations |7T| is eliminated. 
Figure £0] shows a discontinuous behavior of the condensate at T c . Such a jump 
is characteristic of first order phase transitions. This result differs from the one 
in [71], where the transition was found to be second order, although at the same 
T c . Thus, inclusion of thermal fluctuations with self-consistently resummed masses 
renders the phase transition first order. We also show in figure |4.1| how the discon- 
tinuity of the order parameter decreases with decrease of the coupling constant. 
For a weak enough coupling the transition is no longer first order, but a contin- 
uous second-order transition. This can be understood in the following way: the 
self-consistent gap equation ( |4.19|) has been derived in the 0(N) model, which is 
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a generalization of the sigma model for N fields. In such a model, the coupling 
constant should be written as X/N H79fl , assuring the fmiteness of the theory in the 
large- N limit P0fl. We can then write: 



m? = m? — 2— 
N 



3(aj) eq + (N-l)(7rj) eq . (4.22) 



For N — > oo the contribution from the sigma field fluctuation disappears. The 
condensate is then 

- 2 = /. 2 -^(^) eg -/. 2 -^ (4.23) 

In the last equality we have set N—l ~ N for large N. The above expression clearly 
shows a continuous decrease of the order parameter with increasing temperature. 
In the 0(4) model decreasing A by hand is equivalent to going to the large N 
limit in the 0(N) model. We determined a large N critica i ~ 1800 at which the 
transition is second order. For this A = 18 was held fixed, in this way assuring a 
correct vacuum mass for the sigma meson. This result is equivalent to having a 
^critical — 0.01 for which the transition is second order, with N = 4, corresponding 
to the correct number of degrees of freedom. Our results are consistent with the 
qualitative behavior found in [ jS0[] . 

Therefore, in the theory with massless pions there is a first order phase tran- 
sition from a spontaneously broken symmetry phase to a symmetric phase which, 
for small A or equivalently for large N, becomes a second-order phase transition. 
The transition point is characterized by the vanishing of the sigma meson mass. 

4.3.2 Explicitly Broken Symmetry 

When the cx-model Lagrangian includes a symmetry breaking term the behavior 



of the equilibrium condensate is described by equation (|3.25| ): 

\v 3 + A [3(crj) eq + {N- l)(7rJ) e J v - \v 2 v - H = , (4.24) 
and the field fluctuations are given by 

w=2?L dp k^^r (4 - 25) 
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and by 



P 



2 



= dp k^^T ■ (4 - 26) 



The masses in E a = Jp 2 + ml and E w = J p 2 + m 2 , derived in section |3.2.2| are 



ml = 2Xv 2 + — (4.27) 

v 



and 



ml = -= f -ml , (4.28) 

v v 



so that (|4.25|) and (|4.26| ) are the one-loop tadpoles of the dressed sigma and pion. 



Once again, determining the condensate and the meson masses requires a self- 
consistent approach: approximate the value of the condensate v, then solve for 
the masses m a and m n , thus obtaining {cr 2 ) eq and (vr|) e9 , which in turn is used 
to determine v by solving equation ( |4.24| ). Keep repeating these steps with this 



newly obtained v as long as the difference between the old and new values is below 
a required accuracy. The degree of symmetry breaking is given by H = f n m 2 , 
where m n = 138 MeV is the zero temperature pion mass, and v g = f 2 — m 2 /\, 
with A = 18, f n = 93 MeV. With these choice of constants the zero temperature 
sigma meson mass is m a = 575 MeV. 

The numerically determined self-consistent solutions for the condensate and the 
meson masses are displayed in figure |4.2| . In this case a qualitatively different 
behavior is observed. There is no phase transition. The equilibrium condensate 
monotonically decreases with increasing temperature. Even though one cannot 
define a critical temperature, there is a crossover region where the sigma and pion 
mass starts to approach degeneracy. This region, 160 < T < 225 MeV, corresponds 
to much higher temperatures than the critical temperature in the exact chiral limit, 
T c ~ 130 MeV. The reason is that in the present case fluctuations are smaller 
because the masses are larger. The temperature must be at least of the order of 
the pion mass in order to have significant fluctuations which can start destroying 
the condensate. 
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Figure 4.2: Temperature dependence of the meson masses, and of the equilibrium 
condensate in the 0(4) model with explicitly broken symmetry. 
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Chapter 5 

Pion Propagation in Medium 



When particles go through a medium, not only their properties are modified but 
also their dynamics: for example, their velocity is changed. Sigma mesons are 



heavy and their velocity is not altered in any significant manner [|Slq| . Pions, on 
the other hand, massless or massive, have their speed altered by the hard modes. 
In this chapter the change in velocity is determined. First, Goldstone modes are 
analyzed in coordinate space. Then the modification of the pion dispersion relation 
compared to the vacuum relation is discussed, taking into account the resummed 
meson masses. Finally, the dispersion relation of massive pion waves is determined 
for different temperatures. 



5.1 Goldstone Bosons 

In the limit of exact chiral symmetry pions are massless. At zero temperature 
Lorentz invariance requires that massless bosons travel with the speed of light. 
However, at finite temperature, the existence of the thermal medium singles out a 
reference frame, breaking the Lorentz invariance. A direct implication of this is the 
modification of the speed of the Goldstone bosons. This idea has been previously 



looked at by Pisarski and Tytgat in [31]. They studied the propagation of cool 



pions through a thermal medium. Cool means that their discussion is based on 
the low temperature expansion of the self-energy. The analysis presented here is 
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different. We first evaluate the physical response functions directly in coordinate 
space and show that our results are consistent with the results of [pl|. Next, the 
self-energy of the soft pion field is expanded in frequency and momentum to up 
to fourth order. The starting point of our analysis is, once again, the pion field 



equation (|3.17|) : 

d 2 n s + 2\v5(a f n f ) + m\-K s = . (5.1) 

5.1.1 Coordinate space: expansion of the field about its 
equilibrium value 

The effect of the hard modes on the velocity of the soft modes may be seen by 
renormalizing the equation of motion. The idea is to expand the slow field tt s (x') 
within the response function, 5 (a fir f), in a Taylor series about the point x. This 
can be done since we assumed small deviations from equilibrium. The expansion 
is in all four coordinates. The general form of a four- dimensional Taylor expansion 
is 

f(a + x,r + t) = f(a,r)+tf(a,r) + t -f(a,r)+tx-Vf(a,T) + t -x-Vf(a,r) 
+ I(f-V) 2 /(5,r) + ^(f-V) 2 /(a,r) + ... (5.2) 

For now, let us disregard possible dissipation and keep only terms that are relevant 
for renormalization of the equation of motion: 

ir s (t', x ') ~ ir s (t, x) + H s (t, x) + - [{x - x' ) ■ V] 2 7r s (t, x) (5.3) 



The result of inserting the expansion ( |5.3D into the response function ( |4. and 



this into ( |5.1| ) is the renormalized equation of motion: 

tt s (x) — V 2 7r s (x) + m 2 7r s (x) + cii7i s (x) + a 2 Tc s (x) + a 3 V 2 7T s (x) = 

(1 + a 2 )7f s (x) - (1 - a 3 )V 2 7r s (x) + (m 2 + a x )'K s {x) = 
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TT s (x) - V W V 7T s (x) + m 7T 7T s (x) = . 



It has been proven in section [4.2| that the pion mass, which includes one-loop 
corrections, is zero (Goldstone's Theorem), 



-2 mj + ai 

™ v = -T7 = • 

1 + a 2 



(5.4) 



Thus the equation describing the dynamics of Goldstone modes is a familiar wave- 
equation: 



itJx) -vlV 2 TiJx) = 



(5.5) 



where the speed of propagation of the pions is given by 

2 1-03 



1 + a 2 



03 — 02 • 



(5.6) 



In the following we evaluate a 2 and a 3 explicitly^]. 



a 2 



i\m 



,T-,2 



L * <* - *'> 2 /(&/<&/ " )F « * 

2 1 d 3 p 1 (l + ^ + f^E.-E^ + iU-UXE. + E^f 



Am 



A 



ml 



(2n) 3 E (7 E n 
d 3 p 



{El - Elf 



(2tt)' 



i + 2A)^±^-(i + 2/;)^ +: ^ 



(5.7) 



Here E a = yp 2 + m\ and = p. Determining 03 is more cumbersome because it 
requires the extraction and evaluation of the V 2 7r s term from [(x — a?) • V] 2 7r s of 
the Taylor expansion ( |5.3|) . In Cartesian coordinates 



(f- V) 1 



dz 



x 



' d d d 
, ox ay 



dx' 



dy 2 



d 2 



+ z z — + 2xy 



d 2 



dxdy 



+ 2xz 



d 2 



dxdz 



+ 2yz 



d 2 



dydz 



lr rhe lower limit of the time integration is extended to —00. The arguments for doing so are 



presented in section 4.2 
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One can neglect the cross-terms, since these do not contribute to V 2 . The retained 
expression is 



(5.8) 



(2tt) 3 J (2vr) 3 

x(«B- 2 + (y- y + (z- z ™ 



dx 2 dy 2 1 dz 2 

Then 2 separate the integrand as F = j a F a + f^F^ + F , with 

p _ 1 ( e i{E a +E n ){t-f) _ e -i(E a +E„)(t-tf) _ e i{E a -E 1T )(t~t') + e -i{E a -E 1T )(t~t i y 

4E a E n V 

P _ 1 fj(E a +E n )(t-t') _ -i(E„+En)(t-t') i i(E a -E n )(t-t') _ -i(E a -E v )(t-t') 

1 ^(E^+^Xt-f) _ e -i(S CT +E w )(t-t') 



and use the following formulas, when evaluating the integrals: 

dxf(x)5(x - a) = f(a) , (5.9) 



Expression ( |5.8| ) is rewritten as: 

—i\m 2 / (it' 
Jo 



rf 3 p / gVg 2 ^ | d 2 (j)d 2 F (7 | d 2 (j)d 2 F u 
(27r) 3 ct \<9x 2 <9g 2 (9y 2 <9g 2 cfe 2 <9g 2 



(5.10) 



q=-p 



l°° dx> J°° dy' J°° dz ' e ^'+k vV '+^') (>g + ^g + 2 , 2 |V) 



a 2 

9s 2 



| / cfeW fea5 ' / dy'e^' [ dz'e lk * z ' + ^ / da/e****' / dy'y' 2 e ik ^' I dz' ^ 



d 2 (/> 
dy 2 



-(2?r 



|^ y daZe**-*' /" dyV^' / dz'z' 2 e lk ~~ 
^d 2 (j)d 2 8{k x ) 

■ / i i ;.■ km i -4— 

<9y 2 



9a; 2 dk 2 . 



8(k y )S(k z ) + ^ S (k x )^^5(k z ) + ^6(k x )5(k y ) 



dk 2 y 



dz 2 



d 2 6(k z 
dk 2 
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+ 



d 3 q fd^d^ cN^&F^ d 2 (f)d 2 F 7T 



+ 



(27r) 3 V <9x 2 <9p 2 dy 2 dpi dz 2 dp 2 



d 3 p (d 2 <pd 2 F <9 2 0<9 2 F O d 2 ct)d 2 F 



(2ir) 3 \dx 2 dq 2 dy 2 dq 2 dz 2 dq\ 



Notice that Fq, F a and F n are quadratically dependent of the momentum, and so 
these are symmetric in p x , p y and p z . Therefore, the derivatives with respect to 
the different components of the momentum are equal, and one can factor these 
out. The expression of interest is then 



«3 



—iXm^ 



-\m 2 



dt' 



d 3 p 

J a <s n T 



(27T) 



1 



d 3 q cPF^ 

7Jn r% n ~r 



_d^f^^( 

(2nyE a dq 2 x \E»-E» i 

d 3 p 1 d 2 ( '. 



d 3 p d 2 F 

(2vr)3- dp 2 x ' J (2tt)3 dq 2 

_fq_f^f_ ( 1 A 
(2rr^E w dpl \El-El) 



(2nf2E a dq 2 x \E„(E ff + E„) 



(5.11) 



Here E a = yp 2 + m 2 a and E n = q. The last term is the vacuum contribution and 
is divergent for large momenta. So is the term proportional to 1 in ( |5.7| ). These 
divergences are removed with a usual vacuum renormalization procedure. For the 
purpose of our study here, let us assume that this has been done and ignore all 
these terms, focusing on the finite temperature contributions only. 



a 3 



-Xm 

2A 

m 2 



d 3 p fa 

E„ 




d 3 p 
(2^)3 



3 m 2 -EL V 3 m 2 



Putting the ingredients ( |5.7| ) and (|5.12|) into ( |5.6|) one obtains: 



(5.12) 




m 2 
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At low temperatures the large sigma mass leads to exponential suppression of 
the term involving f a = \j[e Ea l T — 1) in (|5.13|) . The solid line in figure |4~T 
shows that the sigma mass is large up to quite high temperatures. Therefore, this 
approximation breaks down only very close to the critical temperature. The rest 
of ( |5.13|) can be done analytically since the pion mass is zero: 



dpp 2 



E n 3 ml 



16 

3ml Jo 



dp 



er — 1 3ml 15 



16 7T 4 T 4 

■>-i 2 



(5.14) 



The final expression for the corrected pion velocity is 



16 



-Xtt 2 



45"" mi 



(5.15) 



Equation fl5.15|) describes the expected result: Goldstone waves travel slower than 
the speed of light when propagating in a thermal medium. Our result is consistent 
with the result in [81|; however, it should be emphasized that our only assumption 
is that the condensate is only slightly out of equilibrium. 



5.1.2 Momentum space: expansion of the self-energy about 
zero frequency and momentum 

Second Order Corrections 

In this section we present an alternative way for obtaining expression (|5.15| ) for the 
velocity of soft pion modes. The propagation properties of particles in a thermal 
medium are governed, in general, by the position of the pole of the propagator. 
This is determined by Fourier transforming the equation of motion ( |3.17|) into 
frequency-momentum space. One readily obtains 

kl = k 2 + m w 2 + U n (k ,k). (5.16) 
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Here (ko, k) is the pion four-momentum. The pion self-energy, U n (ko, k), has a real 
and an imaginary part. Taking k real means that the frequency has an imaginary 
part, ko = u — iT, too. We assume that the imaginary part, T, which determines 
the width, is small relative to the real part, u, which is the energy. Thus, the 
pions are only slightly damped. The question of damping is addressed in detail 
in chapters || and [7|. For now, let us keep our focus on the real part of equation 
( |5.16| ). This is the function relating the energy of the soft pion, u, moving through 
the heat bath, to its momentum, k, known as the pion dispersion relation: 



k 1 + ml + ReU n (u, k) . 



(5.17) 



The pions in the condensate, being soft, allow one to expand the real part of 
the self-energy about zero frequency u = and momentum k = \k \ = 0. The 
dominant non- vanishing terms are 



ReU n (u,k) ~ Re 1^(0,0) + 

d 3 p 



ld 2 Re 11^(0,0) 2 lci 2 Re 1^(0,0) 



doo 2 



2\ml 



(2tt) 3 
bx + b 2 u 2 + b 3 k 2 



G(u,k) 



dk 2 



k 2 



(5.18) 



where 



G(u,k) 



4E a E n 



(1 + fa + U 
+ (U~fa) 



1 



1 



uj — E a — E m u + E a + E, 
1 1 



u — Ei q + E n uj + E a — E w 



(5.19) 



Equations ( |5.18| ) and ( |5.17| ) together lead to 

u 2 (l-b 2 ) = k 2 (l + b 3 )+m n 2 + h 



and so 



OJ 



vlk 2 



ml , 

7T ' 



(5.20) 
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where we defined 



and 



vt ~ 



m 2 = m 2 + 61 . 

7T 7T ' 1 



'3 , 



(5.21) 



(5.22) 



In the following let us determine b\, b 2 and 63 explicitly from their definitions. The 
mass term is 



h = Ren w (0,0) = 2Am 2 |-|^G(0,0) 



2\ml 
2A 



d 3 p 



(2tt)* 



4£L 



(27r) 3 4K^ V El-E, 



-±E a 



El -El 



r d 3 p 


7. 


.A 


1 (2tt) 3 


_ -Eq- 


E w _ 



(5.23) 



The result ( |5.23| ) inserted into ( |5.22j ) together with expression ( |4.9| ) for m^ yields 
another proof of Goldstone's Theorem at one-loop order 







(5.24) 



Relation (|5.20 ) shows the existence of zero frequency excitations at zero momen- 
tum. 

Straightforward calculation leads to 

d 3 p d 2 G{0,0) 



1 (pReUJO, 0) ,_ 2 
2 ^ = 



Am 2 



2A 



mi 



d 3 p 
(2^)3 

<i 3 p 



Vo 



(2tt) 3 



du 2 



E a 



E w 



mi 





1 






(3 + 4- 




E a \ m%) 



(5.25) 



When determining 63 we use the separation ( |5.19| ) 



ld 2 Ren„.(0,0) 



dk 2 



Am, 



d 3 p d 2 G{0,0) 
(2tt) 3 dk 2 



\m„ 



d 3 p 
J2^) 3 



fGA^O) d 2 GV(0,0) 

JO 7; 9 I Jtt 



dk 2 



dk 2 



(5.26) 
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In the first term E a = yp 2 + m 2 a and E w = y p 2 + k 2 — k ■ p, whereas for the second 

term a change of variables has been done, resulting in E a = Jp 2 + m 2 + k 2 + k ■ p 
and E n = p. Then 



h = -ir 



2A r d 3 p 



m 2 



(27T) 



eA 1 + 3 m 2 ) + E^i 1 3 m 2 , 



(5.27) 



Notice that b 2 = —a 2 and b 3 = — a 3 , given by expressions ( |5.7| ) and ( |5.12| ) re- 
spectively, resulting in the same dispersion relation for the Goldstone modes, as 
obtained in configuration space, 

U 16 P 2 



1 



A 



ir 2 m 2 Jo 



dpp 



_Af 4 + ^ 

E n 3 m 2 a E a \ 3 m 2 a 



(5.28) 



The vacuum dispersion relation for massless pions, which are actually on the 
light cone, is 



uj 



(5.29) 



In sections |5.1.1| and |5.1.2| the second order correction to this relation has been 
determined. We call the correction second order, because it is of the order of k 2 . 
The corrected dispersion relation is 



to = v-^k 



(5.30) 



where v n is the solution of 



7T^ 



4 



16 r°° , 
»io dpP E 



m 2 Jo 



dpp 



2 fa 

~eZ 



4 fa 



16 r°° , 

mih dpp e, 



16, 2 T 4 
45 



mi 



(5.31) 



This result coincides with the result obtained in coordinate-space, expression (|5.15|) . 



Figure |4.1| shows that the sigma is heavy, T <C m a , right up to the critical temper- 
ature, which justifies the second line of ( |5.31| ). This takes into account that sigma 
terms are exponentially suppressed. The dispersion relation ( |5.30|) is the same as 
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obtained for sound waves propagating with velocity v n in a non-viscous fluid (see 
for example |39|). Therefore, up to second order the phase velocity, 

v p n h = ^ = v n , (5.32) 



and the group velocity, 



< = (5-33) 



of a Goldstone wave are the same, and are decreasing as the temperature is increas- 
ing. This can be attributed to the fact that the higher the temperature, the more 
are the interactions between the different modes, and these can cause significant 
slowing down. 

Fourth Order Corrections 

One can analyze the dispersion relation further by including higher order cor- 
rections to the velocity. The next-to-leading order non-vanishing terms in the 
self-energy expansion about small energy and momenta are of fourth-order. 

Re IL„(u, k) ~ h + b 2 u 2 + b 3 k 2 + b 4 u 2 k 2 + b 5 k 4 + b 6 oo 4 . (5.34) 

To determine the energy in function of momentum one needs to solve a quadratic 
equation in u 2 , equation that results from inserting ( |5.34| ) in (|5.17|) . The result is 



to 2 ~ k 1 



l + b 2 + b 3 + k 2 {b 4 + b 5 ] 



= vl(k)k 2 . (5.35) 
An analytic result for v v can be obtained when m a 3> T. The result is 



16, 2 T 4 k 2 (XT 2 4Xn 2 T 4 \ 

v 2 w (k) = 1 - — Xjc 2 — + — -— + — j . 5.36 

45 ml ml \ 6 ml 5 mi 1 



When including fourth order terms, the velocity depends not only on temperature, 
but also on momentum. Therefore, the group and phase velocities are different. 
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5.1.3 Numerical Results 



The temperature dependence of the speed of massless pions is presented in figure 
|5.1| , where four different expressions for this are compared. The following notation 
is used: 

v x = 1 (5.37) 

is the speed of light (solid line); v-i is the result of our self-consistent calculation 
given by equation Q5.31J ), using temperature-dependent sigma mass, m 2 = 2Xv(T) 2 , 
determined in Chapter f|, (dashed line) 

vl = 1 - ^Att 2 ^- ; (5.38) 
2 45 fn% y ' 

then 

4 = l - (5.39) 
6 45 ml 

takes into account the temperature-dependence of the coupling Xv(T), but assumes 
the T = mass for sigmas in the heat bath, m a = m a ~ 600 MeV (dotted line); 
finally 

8-7T 2 T 4 

45 j^mi 

is equation (24) from ||81|| , where both the coupling and the sigma mass is set to 
their T = values of v = f n and m a = m a (dotted-dashed line). For all of these 
possibilities one can clearly see a deviation of the pion velocity from the speed of 
light at finite temperatures. 

Let us look at \jfa A a in v% more closely. Here the resummed sigma mass 
m 2 = 2Af 2 introduces not only temperature-dependence but also coupling con- 
stant dependence. For small coupling, A — » 0, v 2 is blown up by \jm\ ~ 1/A. 
Such behavior is not encountered in v 4 . Does this mean that the effect of the 
medium cannot be included in a perturbative manner and a non-perturbative ap- 
proach is needed? One way to reconsider this situation is by looking at v%. Here 
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Figure 5.1: Temperature dependence of the speed of massless pions. See text for 
definition of vi, v 2 , v 3 and t> 4 . 
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the temperature-dependence only of the condensate is taken into account, leading 
to a A 2 correction to the velocity. This does not introduce any bad behavior, since 
with no interactions, A — > 0, the massless pions travel as they should, with the 
speed of light. 

The temperature dependence of the dispersion relation of Goldstone waves 
is presented in figure |5.2| . All the terms of the 2nd-order (|5.31|) (solid curves), 
and the 4th-order ( |5.36| ) (dashed curves) velocity are included in the numerics. 
Also, the resummed sigma mass determined in section is used. This figure 
shows a softening of the dispersion due to second-order corrections compared to 
the dispersion relation in vacuum. This softening becomes significant at higher 
temperatures, and is more accentuated as the temperature increases. The 2nd- 
order curves are straight lines whose slope gives the velocity of the massless pions. 
The deviation of this from the speed of light is clearly noticeable. The 4th-order 
contribution to the velocity is opposite in sign (see (|5.36| )), so its inclusion results in 
a harder dispersion relation compared to the 2nd-order results, but still softer than 
the vacuum dispersion. With increasing momentum a monotonically increasing 
behavior of the 4th-order contribution with respect to the 2nd order one is found. 
The slope of these curves determines the group velocity. 

The behavior of the group velocity in terms of momentum for different temper- 
atures is shown in figure One can see that 4th-order contributions introduce 
a momentum-dependent change in the velocity compared to the 2nd-order results. 
The deviation from the speed of light is greatest for true Goldstone bosons (modes 
with k — 0). Looking at harder excitations, but still in the soft sector, one finds a 
decreasing deviation from the speed of light. Notice that at a certain momentum 
the group velocity becomes super luminous, setting an upper bound for the valid- 
ity of our model. This limiting momentum, denote it by A c , is shifted to smaller 
values for higher temperatures. The scale A c is controlled by the coefficient of 
the k 4 term in ( |5.36|) . The overall validity of the self-energy expansion is limited 
to momenta as large as 125 MeV. This A c can be thought of as the momentum 
scale defined in section |3.2| , which separates hard modes from soft modes. One can 
claim, then, that only modes with | k \< A c can be treated classically, as part of 
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k [MeV] 



Figure 5.2: Pion dispersion curve in exact chiral limit for T = 110 MeV (upper 
curves) and T = 130 MeV (lower curves). 
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Figure 5.3: Group velocity and phase velocity of Goldstone-waves versus momen- 
tum for T = 110 MeV (upper curves) and T = 130 MeV (lower curves). 
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the condensate, and this is consistent with the original premise of our model. 



In figure |5l| the momentum dependence of the phase velocity is also plotted. 
This is related to the index of refraction of the medium, n = c/v, showing that the 
heat bath is less dispersive for Goldstone waves of larger k. 



5.2 Massive Pions 



For the theory with explicitly broken chiral symmetry the dispersion relation, 
determined by Fourier transforming equation ( |3.28| ), is given by 

c 2 = k 2 + ml + ml + Re U n (uj,k) . (5.41) 

Here m^ is the pion mass given by expression ( ET29D in terms of the T = pion 

mass, m n , 



ml 



(5.42) 



and rhn is an apparent contribution which cancels out later. Knowledge of a more 
accurate pion dispersion relation is expected to lead for instance, to more accurate 
description of the dilepton spectrum in heavy ion collisions ||83|| . The method for 
determining an explicit relation between the energy uj and momentum k is the 
same as in the previous section: Taylor expand the self-energy about small energy 
and momentum, retaining the nonzero leading (2nd) order and next to leading 
(4th) order terms, as in fl5.34| ). Because it is lengthy, the evaluation is omitted, 
and only the final expression for the dispersion relation is stated: 



lu = \Jvl(k)k 2 + ml 



(5.43) 



with the function 
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(5.44) 
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where 

h = / fr/«r, h = / |r/a, /a = / |r/., ^ = / fr/* • (5-45) 

JO ti a Jo ti a Jo rb-rr JO Pj^ 



The energies, E a = J p 2 + m 2 and E w = yp 2 + m 2 , are evaluated with temperature- 
dependent masses determined in section |4.3| . One should be careful not to identify 
as a velocity. Due to the nontrivial nature of the dispersion relation, finding the 
group velocity requires the evaluation of du/dk, or looking at the slope of curves. 

The temperature dependence of the dispersion relation of massive pions can be 
obtained only by performing numerical computations. Figure |5.4| shows dispersion 
curves corrected to 2nd-order in k for different temperatures. We found that with 
increasing temperature there is a f lattening of the curves, which corresponds to 
decreasing group velocity. At around T = 180 MeV the velocity of all the excita- 
tions is vanishing. Going to even higher temperatures results in a negative slope 
and so this is a clearly unphysical region. Inclusion of higher order contributions 
changes the qualitative behavior, as shown in figure ETR At temperatures as high 
as 200 MeV there is not only negative slope at low momenta, but also a very steep 
rise at higher momenta, leading to speeds that exceed the speed of light. As before, 
at temperatures this high our results are unphysical. This temperature range is 
the same as the temperature range where the difference between the sigma and 
pion mass is diminishing, and so the pion starts becoming degenerate with the 
sigma (see figure |4.2| ): symmetry is approximately restored. Our result could be 
signalling that at temperatures above 180 MeV pions are not good degrees of free- 
dom anymore. A comparison of results with 2nd- and 4th-order contributions are 
presented in figure [5.6|. These are compared to the vacuum dispersion relation, too. 



Here vacuum dispersion means v n = 1, and sow= \Jk 2 + m 2 . Deviation from the 
vacuum dispersion relation is significant above T = 100 MeV. In the very soft 
sector there is no noticeable difference between the two corrected curves. Their 
behavior is both in qualitative and quantitative agreement with corresponding re- 



sults obtained in a recent publication fl82|. The ground for comparison is limited, 
though, because Ayala et al. studied the kinematic region with k < O^m^ and 
temperatures of the order of the pion mass only. Significant difference between 
the 2nd- and 4th-order curves becomes noticeable in the kinematic region with 
k > 0.3m T . 
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Figure 5.4: Dispersion curves for massive pions corrected to 2nd-order for different 
temperatures. 
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Figure 5.5: Dispersion curves for massive pions corrected to 4th-order for different 
temperatures. 



75 




Figure 5.6: Pion dispersion curve oj 2 = v 2 {k)k 2 + fh 2 for different temperatures 
in the theory with explicitly broken chiral symmetry. Results with v n — 1 (solid), 
v 2 — a (dashed) and v 2 = a + bk 2 (dotted) are shown. 
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Chapter 6 

Dissipation at One-Loop Order 



Dissipation of the condensate occurs because energy can be transferred between 
the condensate and the heat bath through the interactions of slow and fast modes. 
The dissipation can be determined by analyzing how interactions change the width 
of particles or collective modes, since this is related to their damping. In the 
present chapter the effect at one-loop order is calculated by evaluating the response 
functions. We show the direct connection to the physical processes responsible for 
dissipation. After a general discussion, damping of the long-wavelength sigma and 
pion fields is determined, and numerical results are presented. 



6.1 General Formalism 

The equations of motion for long- wavelength sigma and pion fields (eqs. (|3.16| )and 
( P-17| ) in the exact chiral limit, and eqs. Q3.26] ) and ( |3.28| ) for explicitly broken 
symmetry) take the form 

d 2 <p(x) + m 2 0(x) + F(x) = , (6.1) 



where = a s or n s . Based on the response functions (see (|3.42|) , (|3.43| ) and Q3.44Q ) 
one can write 

F(x) = [ d 4 x'(j)(x')n(x - x) . (6.2) 
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Transformed into frequency-momentum space equation ( |6.1| ) reads 

-k 2 + m 2 + U(k) = . (6.3) 
The function Tl(k) is given by 

n(fc) = -if I rfV f^jf 7 ^- C *(*-^«)(*-*')p 1 (p)p 2 ( g )(l + /(pO) + /(g0)) 



<? 2 



(2tt) 4 7 (2tt) 

« 3p 1 r ii + / I+ /y 1 



(2?r) 3 

^ (a; - Ei + E 2 + ie ,c E x - E 2 ic 



uj - Ei - E 2 + ie uj + Ei + E 2 + ie 



(6.4) 

where g is the coupling constant of cubic interactions, and the indices 1 and 2 refer 



to either of the hard modes, <7f, and 7Tj, respectively. Also E\ = \Jm\ + (p + k) 



and Ei = ym 2 + p 2 , and k = (k°, k) is the 4-momentum of the soft sigma meson 
or pion. The frequency has a real and an imaginary part, k° = uj — iT, assuming 
that k is real. The last equality in (|6.4j ) was obtained with the insertion of the free 
spectral functions (see discussion in section |4.2|) . Equation ( |6.4| ) coincides with 
). This is a clear proof that the effect of the medium as described by H(k) can 



be identified with the self-energy. 

The real part of the self-energy participates in the dispersion relation 

uj 2 = k 2 + m 2 + Re n(w, k) , (6.5) 

which has been analyzed in previous chapters. With the assumption of weak 
damping, r uj, the imaginary part of the self-energy completely determines the 
width of excitations: 

- Imn( ^ } . (6.6) 



2uj 

T is the rate at which mesonic modes with energy uj and momentum k approach 
equilibrium. This rate is determined by physical processes which can be identified 
from the imaginary part of the self-energy. 

There are several diagrams contributing to the self-energy at one-loop order. 
The tadpoles are real and their effect is only to modify the mass (see chapter |]). 
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Dissipative effects come from non-local diagrams like the one shown in figure [B.l . 
The imaginary part of this Feynman diagram is given by equation ( B.8 ). There 
are eight processes embedded in ( |B.8| ). For details the reader is referred to [Q. 




<p 2 




01 



a. 





Figure 6.1: Reactions with uj > and mi > m%. (a) requires 
lo 2 — k 2 > (mi + 7712) 2 , and (b) possible for m\ — m\ < lo 2 — k 2 < 
(mi — m^) 2 . 



Figure |6.1| shows the two pairs of diagrams that are possible for positive en- 
ergies, lo > 0. In the following we focus on these processes, since such processes 
are responsible for the appearance and disappearance of pions and sigma mesons 
from the condensate. For uj > the contributions to the imaginary part of the 
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self-energy are limited to 



— * Cy 

Im Il{uj,k) = —7rg 



I 



d 3 p 1 



[(1 + fi + / 2 )% -E 1 - E 2 ) 



(2tt) 3 AE x E 2 



+(f 2 -f 1 )5(uj-E 1 + E 2 )} . (6.7) 



For fixed masses and given values of the energy and momentum only one of the 5- 
function constraints can be satisfied. For the sake of definiteness, assume nix > m 2- 
The process in fig. |6J](a) with u> = E\ + E 2 is possible for u 2 — k 2 > (mi + m 2 ) 2 , 
and the pair in fig. |6.1| (b) with u + E 2 = Ex is possible for m\ — m\ < uj 2 — 
k 2 < (mi - m 2 ) 2 . Rewriting 1 + fi + f 2 = (1 + /i)(l + f 2 ) - fif 2 allows for 
identification of the physical processes directly. These include the decay — > 
4>i4> 2 with statistical weight (1 + /i)(l + f'2) minus its inverse, which is the fusion 
0102 — y ; with statistical weight fif 2 . Note that fusion is an exclusively finite 
temperature process. Similarly, write f 2 — fi = f 2 (l + f\) — f\{l + f 2 ), and identify 
the corresponding reactions as creation 002 — > 0i and decay 0i — > 002- 

In the following we will find explicit expressions for the imaginary part of the 
self-energy and discuss the two contributions separately. 

6.1.1 Processes Above the Two-Particle Threshold 

In the kinematic region where 



forward but tedious. To simplify our task we only consider two cases that are of 
particular interest for our model. 



uj 2 — k 2 > (mi + m 2 ) 2 






Evaluating ( |6.9|) is straight 
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Case I: mi = m 2 = m. 

When evaluating Im II (u;, k) we write the integration measure in the usual way 

/rco rl 
d 3 p = 27cj p 2 dpj dcos6 , (6.10) 

and perform the integral over angles first. This requires a change of variables in 
the delta function. Based on the formula 



(6.11) 



one can write 

5(uj-E 1 -E 2 ) = 5[uj 



dx \ x ~ x i 



2 + {k+p) 2 - yjm 2 + p^j 
5 (uj — \jm 2 + k 2 + p 2 + 2pk cos 6 — ^m 2 + p 2 ^ 



Ei 

= - — (5 (cos 9 — cos 9 ) 
kp 



where 



cos 0n 



uj 2 — k 2 — 2ujy/m 2 + p 2 
2pk 



(6.12) 



(6.13) 



Here k and p are the magnitudes of the 3-momenta. There is a lower limit, A c , for 
the hard loop momentum integration. However, the inequality — 1 < cos#o < 1 
must hold, and this sets natural restrictions on the limits of integration: 

hn U(u,k) = -^J* + dp^-[l + f(uj-E 2 ) + f(E 2 )] (6.14) 

with 



k uj 
P+= 2 + 2 



Am 2 

1 - — and p_ 



uj 2 — k 2 
After a change of variables: 

Im U(cu, k) = - 



k uj 
2~2" 



1- 



Am 2 
uj 2 — k 2 



9' 
16nk j 



r + dE[l + 2f(E)] 



(6.15) 
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where 



■p± + m 2 



2 2 ] 



Am 2 



(6.16) 



cu 2 — k 2 

Provided that uj 2 — k 2 > 4m 2 , where m is the mass of the particle in the loop, one 
gets 

,2 



Im n(a>, k) 



9 



16n 



Am 2 T , 

1 -CT + 2 l log 



1-e" 



1-e" 



(6.17) 



Notice that the T 
contribution is 



and T^O contributions are clearly separated. The vacuum 



Im U(u,k,T = 0) 



167T 



4m 2 
oj 2 — k 2 



(6.18) 



The decay of a heavy boson into two lighter bosons (left diagram in figure |6.1| (a)) 
is possible even at zero temperature. This is textbook material. 

The decay and formation processes are obtained by setting uj 2 — k 2 = M 2 in 
( |6.17| ). Here M denotes the mass of the external particle. 



Im n(M, k) 



16tt 
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1 - exp (— ^ - 
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exp 



UJ I K 
2T "I- 2T 
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It should be emphasized that this expression is different than zero only when the 
threshold condition M 2 > 4m 2 is satisfied; that is, the two-particle production 
threshold is below the pole of the propagator of the boson. Recently, contributions 
to this process from multi-loop self-energy diagrams have been resummed [87 . 
However, inclusion of these is beyond the scope of the present study. 

Since our interest is in the dissipation of soft excitations from the condensate, 
we evaluate ( |6.15| ) for small momenta, k < T 

[(«;+- w_)(l + 2/(a;_))] 



Im U(u, k) 



16irk 
9 2 



16n 



4m 2 
uj 2 — k 2 



(l + 2/(w_)) 



Imn(u;,fc,T = 0) 



1 + 2/ 



2TJ 



(6.19) 
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In case the external boson is much heavier than the bosons inside the loop, these 



can be thought of as massless, m ~ 0, and then ( |6.17|) reduces to 



Im n(c<j, k) 



16n 



T , M-e t 

1 + 2- log 

ft \ 1 — e t 



where 



For fc <T this reads 



ImIT(co>, fc) 



167T 



1 + 2— lo; 



k 10 — k 



(6.20) 



(6.21) 



(6.22) 



Case II: mi = m and = 

In the scenario with one massless and one massive boson inside the loop (see figure 
B.l| ), the evaluation of ( |6.9|) follows the same train of analysis as above, resulting 



in 



JmU{u,k) = --^-r dp[l + f(p) + f(u-p)] 



167rA; 



pi 



with 



u 2 — k 2 — m 2 uj 2 — k 2 — m 2 

Pi = 7 7\ aI1( i V2 - 



(6.23) 



(6.24) 



2(oo + k) ^ 2{u-k) 

The threshold condition that must be satisfied is uj 2 — k 2 > m 2 . The final expression 
is given by 



Im H(uj, k) 



g- 



167T 



1 



m 



T (l-e- P T)(l- 



e t 



+ -rlo, 



u) 2 — k 2 k 



(l — e P t^J (l - < * 



(6.25) 



Once again, a non-vanishing zero temperature contribution can be identified, 



Im U(u,k,T = 0) 



16n 



m 



to 2 — k 2 



(6.26) 



Note that the decay of an on-shell massless boson is kinematically forbidden. 
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6.1.2 Landau damping 

Landau damping is a finite temperature mechanism in which particles are cre- 
ated and/or destroyed through scattering with particles from the heat bath. Such 
processes are displayed in figure |0].b and are described mathematically by 

d 3 p 1 



Im n(co>, k) = —ng 
provided the condition 



(2vr) 3 \E X E 2 



(f 2 - h)5{uj -E,+ E 2 ) 



uj 2 — k 2 < (mi — m 2) 2 
is satisfied. The integrals in ( |6.27[) can be completed resulting in 



(6.27) 



(6.28) 



Im H(uj, k) 



The final expression is 



Im H(uj, k) 



167T k 



, e t — 1 \ I e t 

log — + log 



e t 



e t — 1 



(6.29) 



(6.30) 



where 



p±+m 2 2 , 



k k 2_ u 2 + m 2 i _ m 2 2 (J ^( fc 2 _ UJ 2 + m 2_ m 2 )2 _ ^2^2 _ fc 2) 

P± ~ ± 2 u 2 -k 2 + 2 uj 2 - k 2 



(6.31) 



In the following we discuss the two special cases of interest. 



Case I: mi = 1112 = m 

This case requires the 4-momentum of the external particle to be space-like, that 
is uj 2 — k 2 < 0. Expression ( |6.30| ) is simplified to 



q 2 T I- e~r 

Im U(u, k ) = -f- T log 

8% k \l — e — ~ 



(6.32) 
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where 

uj k / 4m 2 

Expression ( |6.32| ) vanishes at T = 0, just as expected. This process is forbidden 
on mass-shell because it would yield tachyons. 



Case II: ni! = m and m 2 = 

In this case one of the particles in the loop is massless and the other has a kinematic 
constraint, uo 2 — k 2 < m 2 . 

Imn( W) fc) = --^-J 00 dp[f(p)-f(u J + p)} 



9 



16irk 



Vl 

dpf(p) 

Vl 

pi 



g 2 i er — e t 



Tlog — s — , (6.34) 



16vrA; I e r -1 



where 



k 2 — lo 2 + m 2 m 2 

Pl= 2(u + k) ^Jk aS U ^ k - (6 - 35) 

For the self-energy of massless on-shell bosons the result is 

a 2 rinr+ k 

lmU{k)= -4^J^ dpfip) - (6 - 36) 



4k 

When the massless on-shell bosons are soft, too, then 

6.2 Sigma Decay 

Contributions to the imaginary part of the sigma meson self-energy come from 



two different couplings. These were derived in section |3.3.1| , and are 3Xva s aj and 
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Xva s iT 2 . The corresponding coupling constants are g = 3Xv and Xv, respectively. 
The corresponding diagrams are presented in figure |6.2| . The factor 2 in front of 
the diagrams is a symmetry factor due to the interchangeability of the lines inside 
the loop. 




Figure 6.2: One-loop self-energy contribution to the long- wavelength 
sigma a s with couplings to the hard <7/ and -Kf through 3Xva s cr 2 and 
Xva s 7r 2 . 



Let us look at the on-shell processes that are allowed since these determine the 



dynamics of the decay. Figure |J72] shows that mi = rri2 = m, where m = m a or 
m = m n . For an on-shell sigma meson, u 2 = k 2 + m 2 , condition ( J6.28[ ) is never 
satisfied. Moreover, for m = m a condition ( p.8| ) is not fulfilled either. One can 
say that the kinematics forbids Landau-damping and a s — > er/er/ decay. However, 
contribution from the decay of soft sigma mesons into hard thermal pions is non- 
zero above the threshold m 2 > 4m 2 . The rate for a s — > 71771-/ is given by equation 



R9) 



where E\ = ym 2 + p 2 = p and Ei = ym 2 + (p — k) 2 . To obtain the total rate 
one has to take into account the reverse process, the annihilation of two pions from 
the heat bath into a low momentum sigma, 717717 — > a s , which has a rate ( p.!0| ) 



sc. 



The total rate is then 



F(u,k) = T d (u,k) -Ti(uj,k) 
d 3 p 



2. J (ar). AETe^ 2 ^ ~ £l " + h + «■ (6 - 38) 

The scattering amplitude for this process is basically the coupling constant asso- 
ciated with the vertex, \Ai\ 2 = g 2 = (At>) 2 . Performing a simple variable change 
p — k = q in (|6.9| ), the identity 

T(u,k) - 



Im H(u, k) 



(6.39) 



becomes transparent. It is important to be aware of the difference of a factor of 2 
between expressions (|6.6| ) and ( |6.39| ). The first one describes the rate of decay of 
the amplitude of the wave, exp(— Tt/2), while the second one represents the loss 
rate for the number density, exp(— Ft). 

The imaginary part of the self-energy of a soft, fc < T, on-shell sigma meson 
is given by eq. ( |6.19|) as 



Imn» = Imn ff (w 1 T = 0) 1 + 2/ 
where the non-zero vacuum contribution is, from ( |6.1S| ), 

X 2 v 2 



2TJ 



Im U a (uJ } T = 0) 



8tt \ 



Am 2 



mi 



(6.40) 



(6.41) 



The rate at which sigmas of energy u disappear from the condensate due to their 
decay into pions is 



\ 



Mcoth 



-) 

AT) 



(6.42) 



m 



- - 2Xv 2 and that there are — 1 



r ( ,\ - - 1) , m 2 - mj 

167T UJ 

where account is taken for the relation m 2 
pion fields. In the rest-frame of the sigma to = m a . When the pions are Goldstone 
bosons expression (|6.42|) is reduced to 

(N-l) 



[m n 



16tt 



-\m„ coth 



AT 



(6.43) 
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This agrees with equation (23) in |7C|. At the critical temperature m a = and 
just below this m a <C T . Thus one can perform a high temperature expansion 
using 

! oc 

coth(x) = — I h 0(x 4 ) as x -> . 

v ' x 3 45 

When keeping the leading term only, eq. ( |6.43|) results in the classical limit 

r_ ~ ^F^ AT • ( 6 - 44 ) 

1D7T 

At low temperatures the sigma is heavy, T <C rn a , and no such expansion is 
reasonable. In this case 



and 



coth(a;) — > 1 as x — > oo 



r^^(iV-l)^ , (6.45) 



16vr 

which is equal to the decay rate calculated at zero temperature. 



6.3 Pion Damping 

At one-loop order there is only one diagram contributing to the pion self-energy, 
the origin of which lies in the 2Av7r s cr/7rj coupling with g = 2Xv. This diagram is 
shown in figure p.3| . 




Figure 6.3: One-loop self-energy contribution to tt s 
with coupling to a s and n s through 2\v7r s afiTf. 



At this level, dissipation of the pion condensate occurs provided that the energy and 



8N 



momentum of the soft pion satisfies the condition uj 2 — k 2 < [m a — m n ) 2 . Then 
the transformation of a pion into a sigma when propagating through a thermal 
medium can happen. Basically, a soft pion from the condensate annihilates with 
a hard thermal pion producing a hard thermal sigma meson, vr s 7Tj — ► cry. The 
inverse process is the decay of a hard thermal sigma meson into a soft and a hard 



pion, oj — > TTg-Rf. The net rate of dissipation is given by eqs. (B.9|) and (|B.10|) as 



T(u,k) = T d (u),k) -Ti(u,k) 



2uj J (2?r) 3 *E a E v 



\M\ 2 (27i)S(cu + E w -E CT )(f n ~f a ) , (6.46) 



where E a = y m 2 + (k + p) 2 , E w = p, and the amplitude is \A4\ 2 = g 2 = (2Xv) 2 . 
Note that with QrliBD and QQ7| ) 

— * 

Im II(u;, fc) 



r(w,ife) 



The rate of dissipation of massive pions is 

r„,M) = --iA r( '";-'" ;) [log 

o7r /ecu 

where 



e t — 1 
e~ — 1 



+ log 



e t — 1 



e t 



(6.47) 



(6.48) 
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k k 2 - to 2 
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4m 2 . (cu 2 



fc 2 ) 



fc 2 



k m 2 — 2m 2 oj ml 
±- „ - + -- 



4m 2 



cij 2 — A; 2 
for cj 2 = A; 2 + m 2 



2 m 2 2m^ 
The damping rate of massless pions is given by equation fl6.37f ) 



(6.49) 



r (u, k) — — 



1 2rp EL 

-A- — log — 



8tt kuj 



n 
e t 
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where 



Pi 



k 2 — uj 2 + m 2 
2(cu + k) 



(6.50) 



(6.51) 
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On-shell, u = k Goldstone modes are damped as 

Xmt „ I m 2 ^ 



In the low temperature limit T <C m a the approximation 

/(*) = 7 1 



e x - 1 
can be used, leading to 

V ma {k) ~ ^ e -^ 4fcT . (6.53) 

This form clearly shows that the process is exponentially suppressed by the heavy 
sigma meson. As the momentum vanishes, k — > 0, the damping rate goes to zero as 
it should in order to satisfy Goldstone's Theorem. For high temperatures T 3> m a 
and T — > T c 

fix) = ► as x^U 

J w e x -1 x 2 12 



resulting in 



, 1S Xw? (AkT 1\ XT , . 

'W^ssUe-i "aF- (6 ' 54) 



The last expression corresponds to the classical limit. As pointed out in p9fl , 
care has to be taken when taking this limit. For low momentum, k m^/AT, 
the classical approximation would overestimate the damping rate, since this is 
exponentially suppressed. Moreover, in the classical limit the damping is not zero 
when the momentum is zero, thus Goldstone's Theorem is not fulfilled. 



6.4 Numerical Results 

The temperature dependence of the sigma damping rate, eq. ( |6.42j ), in the rest 
frame of the sigma is shown in figure |6.4| . Note that even at zero temperature 
there is a finite damping, which means that the sigma meson can decay into two 
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pions in the vacuum. Our results show that the dissipation is greater in the theory 
with exact chiral symmetry (upper curves) compared to the explicitly broken sym- 
metric case (lower curves). This is due to the fact that there is more phase space 
available when pions are massless than when they are massive. At T = 0, for model 



parameters defined in section |4.3| , the damping rate is about T a7T7T = 600 MeV when 
= and T an7T = 510 MeV when m n ^ 0. These values are of the order of the 
mass of the sigma, meaning that the width of the sigma resonance is very broad, 
in other words the sigma meson is overdamped. Figure |6.4j presents the compar- 
ison of results obtained when calculations were done keeping the masses at their 
zero temperature value (solid lines), and when the temperature-dependence of the 
masses is taken into account (dashed lines). Both set of curves show a suppression 
of the sigma decay rate when the resummed masses are used. This suppression 
is increasing with T. When m n = the decay rate diverges at the critical tem- 
perature on account of large fluctuations. However, this would require further 
investigation. Without this, one can only say that due to its strong damping the 
sigma field quickly relaxes to its equilibrium value. In the ^ case a significant 
qualitative change is exhibited: the increase of the damping with T is followed by 
its quite drastic decrease starting from about T = 150 MeV. This corresponds to 
the temperature where the sigma mass begins to drop significantly (see figure f4.2| ). 
It is then natural to expect a decrease of its decay rate into pions with masses ap- 
proaching that of the sigma. Thus, the strong relaxation of sigma mesons at higher 
temperatures is not obvious anymore and some oscillations in this field might be 
expected. As the temperature drops the damping of the sigma mesons will set in. 



Figure |6.5| shows the temperature dependence of the damping rate of massive 
pions at one-loop order calculated in the pion's rest frame. Unlike the sigma 
meson, at zero temperature the dissipation is zero. This makes sense because the 
transformation of pions into sigmas is due to their annihilation with a hard thermal 
pion in the medium, and so this is an exclusively finite temperature process. Notice 
that T n7Ta <C T a7Tn . This is the case because the phase space available for the pion 
decay process is suppressed by the large sigma mass. As the temperature increases, 
though, the dissipation rate is increasing too. As shown in figure B3I the damping 
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Figure 6.4: Temperature dependence of the sigma damping rate at one-loop order 
calculated in the sigma rest frame in the case of m n ^ (lower) and in the chiral 
limit rrij, = (upper) with T = masses (solid) and T-dependent masses (dashed). 
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Figure 6.5: Temperature dependence of the pion damping rate at one-loop order 
calculated in the pion rest frame with T = masses (solid) and T-dependent 
masses (dashed) and compared to the pion mass. 
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can get quite strong. At T = 170 MeV, for example, when the pion mass is kept 
at its vacuum value, m n = 138 MeV, T n7Ta = 14.4 MeV, and for = 158 MeV 
the damping is T n7Ta = 87.0 MeV. Thus the width of pion excitations is about 10% 
and 55% of their energy, respectively. This result should lead us to rethink our 
initial assumption of T <C u (see section |6~T|) . 



The damping in terms of the momenta of the pion fields is presented in figure |T| 
for different temperatures. Results using the dispersion relation u 2 = k 2 + m 2 with 
zero temperature masses (solid lines) and with resummed masses (dashed lines) 
are displayed. Noticeable damping occurs above T = 100 MeV, and this increases 
with T. One can see that pions that have the vacuum mass and hard momenta are 
more damped than pions with the same mass and soft momenta. However, when 
the temperature-dependence of the mass is taken into account such a statement 
cannot be made. At T = 160 MeV all modes are equally damped. In other 
words, the width of pion modes is independent of their momentum. This width is 
increasing with T and it can be as great as 30% of the energy. We also analyzed the 
momentum dependence of the energy and width of pions with different dispersion 
relations. Our results for the vacuum (solid lines), the 2nd-order (dashed lines) and 
4th-order (dotted lines) corrected dispersion relations (see section |5[) are presented 
in figure [BT^ . No major difference exists between the vacuum and 4th-order results. 
Note, though, that the quantities are normalized to the vacuum pion mass, so the 
true differences are somewhat greater than indicated on the axis. The overall 
message is that the width of the excitations, and so the damping of the fields, is 
more pronounced at high temperatures and should not be neglected. 

We close this section by presenting the results of our study of the dissipation of 
soft Goldstone modes which happens through the absorption of a high frequency 
Golstone boson producing a sigma meson. The momentum dependence of the 
energy and width of these massless excitations at different temperatures ( |6.48 



is 



shown in figure |6.8j . Notice that modes with k = have T = 0, meaning that true 
Goldstone modes do not decay. Thus homogeneous condensates are undamped 
even at high temperatures. This result is consistent with pS|. Figure |6.8| shows 
that the pion modes with higher momenta are significantly damped near T c . 
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Figure 6.6: Momentum dependence of the pion energy (upper) and width (lower) 
at different temperatures calculated with T = masses (solid) and T-dependent 
masses (dashed) and dispersion relation uo 2 = k 2 + m 2 . 
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Figure 6.7: Momentum dependence of the pion energy and width at different tem- 
peratures calculated with temperature-dependent masses and dispersion relation 
oo 2 = vk 2 + m 2 with v — 1 (solid), v — a (dashed) and v — a + bk 2 (dotted). 
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T=110MeV 




Figure 6.8: Momentum dependence of the energy (upper curves) and width (lower 
curves) of Goldstone modes (m n = 0) at different temperatures calculated with 
T-dependent masses and dispersion relation uj 2 = vk 2 with v — 1 (solid), v — a 
(dashed) and v — a + bk 2 (dotted). 
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Our results are contrary to previous beliefs, according to which sigma mesons 
are quickly damped and, due to their small mass, pion oscillations persist in the 
small momentum region [58|. We found that at high temperatures not only is the 
damping of sigma mesons significant, but also that of the pions. This damping is 
most accentuated in the phase transition region. Accordingly, we do not expect 
amplified pion oscillations but rather the decay of non-equilibrium pion fields. 
Clearly, further investigations are required to determine if there are oscillations of 
damped waves surviving. Moreover, to determine how fast thermal equilibrium is 
reestablished the time-scales that are involved should be studied. 
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Chapter 7 

Dissipation at Two-Loop Order 



This section is dedicated to a detailed discussion of processes in which a meson 
from the thermal bath (pion or sigma) elastically scatters off a meson from the 
condensate (pion or sigma). As a consequence of these scatterings mesons get 
knocked out of the condensate, making this analysis relevant for the study of 
dissipation. Contrary to one-loop level damping processes which are present only 
when certain conditions are satisfied, scattering always occurs. First, I present 
the general formalism for determining the rate of dissipation. Then, all possible 
tree-level scattering processes are evaluated. Finally, my numerical results are 
presented. 

7.1 General Formalism 

Two-loop linear response functions can be evaluated following the procedure pre- 
sented in the previous chapter for one-loop order. However, rather than attempt 
this sophisticated evaluation, a direct determination from the physical processes 
is presented. Two-loop contributions correspond to two-particle scatterings with 
amplitudes evaluated at tree level. The equivalence of these two methods at two- 
loop order can be shown in a way similar to the proof at one-loop order presented 
in Appendix [B| 

The contribution from scattering processes to the rate at which the condensate 
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fields decay to their equilibrium values is given in terms of the imaginary parts of 
their self-energies 



TV U\ IlD n ( W ' fc ) ( r 7 1 \ 

— * 

Here u and /c are the energy and momentum of the meson. This rate refers to the 
field amplitude, which is half of what the decrease in number density is. The general 
form of the self-energy of a particle of mass m a , propagating with four-momentum 
k = (u, k) through a medium in thermal equilibrium, is given by dC] 

Here A4 is the transition amplitude for the process ah — > ah. The thermodynamical 
weight f(E) is the Bose distribution of thermal mesons of mass m& and four- 
momentum p = (E,p). In terms of the forward scattering amplitude A4(s) = 
— 87r v /i/ cm (s), where s = (p + k) 2 . From (|7.2|) follows the imaginary part 

Imn tt6 (fc) = - J 0-f(E)^S^a total (s) . (7.3) 

To obtain this equality we applied the standard form of the optical theorem that 
relates the imaginary part of the forward scattering amplitude and the total cross- 
section |73| : 

Im f cm {s) = -^Vtotalis) . (7.4) 

The differential cross section for a(k) + b(p) — > a(k') +b(p') is given by the general 
formula 1731: 



da ab = 1 \M\ 2 dr 2 , (7.5) 



{p ■ k) 2 — rnlm 2 
where dr 2 is the two-body phase-space 



^ = <^ + *-'- V%?w - (7 ' 6) 
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One can simplify this expression by evaluating the phase-space integral in the 
center-of-mass frame of the two colliding particles. Then the total scattering cross- 
section is given by 

\ / cm 

where 

l ' h ' A ' \M\ 2 . (7.8) 



dVL ) 64vr 2 .< 



The symmetry factor 1/51 is due to the number S of identical particles in the final 
state. It is clear that knowing the amplitude j\A of a process readily results in the 
dissipation rate due to that process. However, we must distinguish between scat- 
terings that involve massive and massless mesons. For the former it is convenient 
to work in the rest frame of the meson a from the condensate. Then k — 0, and 
the dispersion is u — m a . It makes no sense, though, to talk about the rest frame 
for massless particles. In this case we evaluate the self-energy in the rest frame of 
the condensate. Then k labels the momentum of the meson with respect to this 
frame. Expression (|7.3|) then takes simplified forms: 

Im U ab (u = m a , k = 0) = / dE[E 2 - m 2 )f(E)a ab (E) , (7.9) 

for massive mesons, and 

lmU ab (u,k) = -— dcosO dp^-(s-m 2 b )f(E)a ab (E), (7.10) 

07T Z J-l JO E 

for Goldstone modes (m a = 0). Here s = ml + IElo — 2pkcos6 where 6 is the 
angle between k and p. The dispersion relation of the hard thermal modes is 



E = \jp 2 + mf. In both of the above formulas there should be a Q(p — A c ) 
restricting the integration to hard momenta only. The probability of finding long 
wavelength modes in the heath bath is small leading to negligible contribution 
from the soft end of the spectrum. 
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7.2 Sigma Scattering 



There are two mechanisms that contribute to the removal or addition of a sigma 
meson to the condensate: elastic scattering of a hard thermal sigma and a hard 
thermal pion off a sigma from the condensate. 



7.2.1 Sigma-Sigma Scattering 

To first order in the coupling A there are four diagrams contributing to the process 
in which a thermal sigma meson knockes out a low momentum sigma from the 



condensate. These are presented in figure [7J]. The transition amplitude is the 
sum of contributions from different diagrams: 



M 



-6A 



+ 



+ 



(7.11) 



s — m£ t — u — m 2 a 

which reflects the symmetry in the s = (pi +P2) 2 , t = (p± — qi) 2 , u = (pi — q?) 2 
channels. These Mandelstam variables satisfy s + t + u = 4m 2 . The total cross- 
section is 



9A^ 

87TS 



s + 2m 2 — 3m 2 



+ 



18(m 2 



m 



2\2 



m^(s — 3ml) 



12(m 2 — m 2 )(s 2 — 3sm 2 — mi + 3m 2 ml) 



In 



3m 2 



mt 



(s — 4m 2 r )(s — 2ml)(s — m 2 ) 
In the low energy limit expand ( |7.12| ) about s = 4m 2 . To leading order 



9 x2 {4ml-hmlf 



32tt 



rrvi 



and in the rest-frame of the sigma this gives rise to 



Im IL 



327T 3 



2\2 



m°„ 



(7.12) 



(7.13) 



(7.14) 



In the high energy limit only the four-point vertex contributes to the amplitude, 
resulting in 

9A 2 



8tts 



(7.15) 
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Figure 7.1: oo scattering at tree-level. Straight lines 
are as, dashed lines are 7rs. 



This gives rise to 



Im IL 



X 2 T 2 . 



64tt 



(7.16) 



With these two limits one can construct an interpolating formula that describes 
the whole energy range: 

i m n g 9A2 ^K + T)(K-5m^ 

327r6(m CT + T)(4m2 -5m2) 2 + 7r 2 m5(e m -/ T -l) ' ( ' J 

In the theory with exact chiral symmetry m w = and (|7.17|) is simplified to the 
result obtained in 

9A 2 T 2 



Im IU l_ 2?r g6 + _ x y 

The contribution to the rate of decay of the amplitude is then 

ImhU 9A 2 T 2 



2m a 4tt m a [96 + vr 2 (e m -/ T - 1)] 



(7.18) 



(7.19) 
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7.2.2 Sigma-Pion Scattering 

A hard pion in the heat bath can be energetic enough to knock out a sigma-meson 
from the condensate. The possible tree-level processes may happen according to 



the four different diagrams shown in figure [7.2| . The transition amplitude obtained 
from these is 



M 



-2A 



l + (m 2 - ml\ 



s — mi 



t — mt 



+ 



u — mi 



where s = (pi + q\) 2 , t = (pi — P2) 2 , u = (pi — q^) 2 and s + t + u 
The total scattering cross-section is 



2(ml 



(7.20) 
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(2m 2 + m 2 — s)((m 2 — m 2 ) 2 — sm 2 ) 
6s(m 2 — m 2 ){sm 2 + 2sm 2 — s 2 + ml — ml — 2m 2 m 2 ) 
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(7.21) 



The low energy limit is obtained by expanding ([7.21 ) about s = (m a + m n ) 2 and 
is 



= 9 x2 mj(3ml - Ami) 2 
47r sm 2 (m 2 — 4m^) 2 



(7.22) 



At high temperatures, where only the first diagram of figure [T2] contributes, (|7.21|) 
reduces to 



Airs 



(7.23) 
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Figure 7.2: an scattering at tree-level. Solid lines are as, dashed 
lines are ns. 



The contribution to the imaginary part of the self-energy is readily determined in 
these two limits using ( |7.9| ). At low energy 

Im IU = -^X 2 T 2 e- m ^ T F an (m a , m n ) , (7.24) 



where 



r cnryfllu, m n ) — „ . ... „ „.„ , [I .ZD) 



and at high energy 



A 2 T 2 

Im IU = -7^7— • (7.26) 
yo7r 



The two limits can be combined into one approximate expression which is used 
then to determine the rate of dissipation 

Q T 2 F 
T ~ —A 2 - — (7 27) 

CT ""87r m CT 216F (T7r + 7r 2 (e^/ r - 1) { } 
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This expression gives the contribution from a single pion. In order to factor in all 



pions, we must multiply (|7.27| ) by N — 1. 

At this point a distinction should be made between scatterings with massless 



and massive pions. Notice that the low energy expression ( 7.22 ) vanishes for zero 
pion mass. The first nonzero term in the series expansion of the cross section is 
the 4th order term, which yields 

,2^4 

(7.28) 



7 



On 



3tt 



-A 



m„ 



giving rise to 



Im IL 



1344 VC(7)11 . 



7T° 



mi 



(7.29) 



The high temperature limit is given by (|7.26|) . We combine the two limits again 
and obtain 



r ~ 



(N 



A 2 T 2 



1927T m a 



T 5 + (rrir/lO.Gy 



(7.30) 



formula identical to the one obtained in WO 



7.3 Pion Scattering 

Dissipation of the condensate can arise from scattering of soft pions with hard pions 
or hard sigma mesons. The damping of massive pions and that of the Goldstone 
pions is expected to be different, requiring a somewhat different analysis. 



7.3.1 Pion-Sigma Scattering 

The possible tree-level diagrams representing the reaction in which a hard thermal 
sigma knocks out a low momentum pion from the condensate are shown in figure 
[7.2| . The total cross section is the same as for sigma-pion scattering and is given 
by expression ( [7.211 ). The low and high energy limits are also the same and are 
given by equations ( [7.22] ) and ( |7.23|) . 
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When the pions are massive the imaginary part of the self-energy in the low 
energy limit in the rest frame of the pion is 



where 



Im IU = - A A 2 T 2 e -m CT /T F7rCT j 



ml(m a + T)(3m 2 - Ami) 2 
m^(m a + m 7r ) 2 (m 2 . — 4m 2 ) 2 



(7.31) 



(7.32) 



and for high energies is 



X 2 T 2 
' 96tt 



Im rL 



(7.33) 



The total scattering rate due to this process can be parametrized by an interpo- 
lating formula between the two known limits: 



r ~ 

- 1 - 7TCT" 



9A 2 T 2 



8nm a 21QF W(T + n 2 (e m °/ T - 1) 
When the pions are Goldstone bosons the cross-section ( 7.21 ) reduces to 



(7.34) 



4tcs 



[s + m 



2\2 



+ 



9sm 2 



2m 2 — s s 2 — sm 2 + vrA 



2sm 2 (s 2 — 3sm 2 — m 4 ,) ^ ( s{2m 2 rT — s) ' 



s(s — m 2 ) 3 



6sm 2 (s 2 — sm 2 — m 2 ) ^ ( s 2 — sm 2 + m\ 



sis — m 



2\3 



smt 



(7.35) 



with low energy limit obtained from the series expansion of this about s = m 2 



7\ 2 (s - m 2 ) 4 



C77 



37r m 40 

When T — > T c the three-point vertex diagrams do not matter and 

A 2 



(7.36) 



Atcs 



(7.37) 
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The self-energies in these limits are obtained from ( |7. 10| ) with s = ml + 2Eu — 
2\p\\k\ cos8. It is perhaps worth noting that when evaluating ( 7.10|) we keep the 
energy uj and momentum k independent. At high temperatures 

A 2 T 2 



Im U ncT (u}, k) 



96tt 



(7.38) 



This is exactly the same as for the massive pion case. No surprise here because 
this limit is defined by the temperature being much larger than the masses. The 
evaluation of the self-energy in the low temperature limit results in 

56v^FA 2 



Im U na (u, k) = 
with the following definition: 

u(uj, k) 



3V2 



'—e- m "/ T u(u,k)uj , 
m„ 



(7.39) 
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lhk A — + lOJfeV 
mi 



- T 2 
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m 2 a 


ml 






ml 



(7.40) 



An interpolating formula for the contribution to the massless pion damping rate 
is: 

A 2 T 2 u(u,k)< 



1927T 



(7.41) 



7.3.2 Pion-Pion Scattering 

7T7T scattering has been extensively studied during the last couple of decades in a 
variety of different models and approaches. An incomplete but significant list of 
references is [|9~2|j-[|9~6|j. 

Here we study the scattering of a hard pion off a soft pion at tree-level. In our 
model with one pion field there are four diagrams contributing to this process: one 
4-point vertex diagram and three 3-point vertex contributions involving a sigma 
exchange in the s, t and u channels. These are presented in figure 7^3 . The 
transition amplitude is 

OS ( 1 1 1 



M 



-2A 



3 + (mi 



+ 



t 



u 



m 



(7.42) 
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Figure 7.3: 7T7T — ► 7T7r scattering at tree-level in 0(2) model. 
Solid lines represent sigmas and dashed lines represent pions. 



where s = (pi + P2) 2 , £ = (Pi — 9i) 2 and w = (pi — 92) 2 , and s + t + u = 4m\. 
Remember that for pion-sigma scattering all pion species were taken into account 
simply multiplying the result by an overall factor of N — 1. In the case of pion-pion 
scattering the situation is more complicated since additional channels open up. Let 
us look at the model with N = 4 for physical pions. With the usual definition of 

7Tl ± 7T 2 



7T 



V2 



and 



7T 



713 



(7.43) 



figure [7l| shows the possible diagrams for 7r°7r° — > -k°7i° scattering. Since we 
are interested in the dissipation of the scalar condensate the cross-sections for 
7r + 7r~, 7r°7r + — > 7r°7r + and 7r°7r~ — > 7r°7r~ processes are evaluated. The 



7T 7T 



possible diagrams are shown in figure |7.4j . The corresponding transition amplitudes 
are 



M 



00^+- 



-2A 



1 + 



ml 



ml 



(7.44) 
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Figure 7.4: 7r°7r° — > 7r + 7r~, 7r°7r + — > 7r°7r + and 7r°7r _ — ► 7r°7r~ scattering at 
tree-level. Solid lines represent sigmas and dashed lines represent pions. The 
4-momenta are shown in brackets. 
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m a — ml 
t — ml 



(7.45) 



Instead of writing long formulas for the cross-sections let us display only the leading 
term of the expansion of these about s = Am 2 . This low energy limit is quite 
acceptable for m a ^> T. The sum of the contributions from all the processes of 
is 



figures [O and 7.4 



A 2 



mt 



. tti 2 Tn^ \ 

, , ^ 23-16^ + 128-J . (7.46) 
327T (m? a - Ami) 2 \ m l m t J 

It is important to mention that this is not the isospin averaged cross-section found 



in other papers p7| . We study only scatterings that involve the neutral pion 
field, because condensates of charged pions would be destroyed by electromagnetic 
repulsion. Therefore contributions from processes like 7r + 7r~ — ► 7r + 7r~ have not 
been taken into account. In the T —> T c limit only 4-point vertices contribute to 
the amplitude, resulting in 

15A 2 



57TS 



In the rest frame of the 7r° the above two limits give 



Im IL 



23 



and 



32tt 3 



Im IL 



A 2 TV 



v/T ml(m n + T) 
4m 2 ) 2 



(mt 



5A 2 T 2 
64vr 



(7.47) 



(7.48) 



(7.49) 



respectively. The rate of dissipation due to massive pion-pion scattering is then 
23A 2 T 2 m 2 (m^ + T) 



r ~ 

7T7T 



64tt 



_5_ 
46' 



Ami' 



(7.50) 



-m%(m.n + T) + 7r 2 (m 2 

When the pion mass is zero the expression ( 7.46Q vanishes. The first nonzero 
terms in the low-energy expansion of the cross-section are 

5A 2 



A 2 s 2 



23A 2 s 3 



12n mi Atc m % „ 107T m % „ 



(7.51) 
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The leading term coincides with the current algebra representation for the total 



cross-section, eqn. (7.14) with m 7T = in [95]. In the rest frame of the condensate 
s = 2Euj — 2pk cos 9 and hard modes are considered to be on the light-cone, E = p. 
Then ( |7.10| ) gives to leading order 

l m U nw (u,k)^-^-X 2 ^(h 2 + u 2 ) . (7.52) 
3o mt \i J 



Expressions (|7.52 ) and ( 7.49 ) can be combined into an interpolating formula. How- 



ever, since the sigma is heavy almost up to the critical temperature it is feasible 
to write the damping rate of Goldstone modes as 

, Im U w Ju, k) 7T , 9 T 4 (Ik 2 \ 
T 7TW (u,k) = P ; ~ —A — h uj . (7.53) 



Expressing dispersion relations in the form uj = v n (k)k, this damping rate vanishes 
linearly as k —>■ in the chiral limit, showing that Goldstone's Theorem is obeyed. 

7.4 Numerical Results 

All the numerics have been done for the model with N = 4. 
7.4.1 m n ^0 

The temperature dependence of the total scattering rate of the sigma meson, 
r aa -\-T a7T , is shown in figure [7~5l . Here expressions (|7.17|) and (|7.27j ) have been used. 



The masses involved are the self-consistently determined meson masses. One can 
see that scattering is more accentuated at higher temperatures. For comparison 
we plotted the energy, which is the sigma mass in its rest-frame. The scattering 
contribution to the sigma damping rates is well below the energy. Moreover, dissi- 
pation of sigmas from the condensate due to their scattering is much smaller than 
due to their decay. This means that the scalar order parameter relaxes to its equi- 
librium value via the production of lighter pion fields. It also means that sigma 
mesons are so unstable that they are more likely to decay before they could ever 
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Figure 7.5: Scattering contribution to the width of the sigma meson as function of 
temperature in the model with explicitly broken chiral symmetry. 
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scatter with other particles from the medium. The total damping rate is obtained 
by adding all the components 



l-V r<T7T7r -|- ^ era 1~ <77T • C^ - *-*^) 

This is dominated by Y U1T1T . 

The scattering of long-wavelength massive pions according to formulas ( J7.34 ) 



and ( 7.50|) is presented in figure |7]6|. Due to the heavy sigma exchange there is a 



strong suppression at low temperatures. When reaching T ~ 130 MeV the scat- 
tering rate becomes significant and it is comparable in magnitude to the damping 
at 1-loop order. 

7.4.2 111^ = 

Damping of heavy sigma mesons due to their scattering with other heavy sigmas 
and Goldstone bosons is shown as functions of temperature in figure [7.7| . This rate 
is more important at higher temperatures, although it is not dominant over the 
damping due to sigma decay. 



In figure [TJI massless pion scattering processes' contribution to the damping 
of the condensate is presented. Notice that all curves start from the origin as 
they should. This figure is the result of expression ( |7.53|) using different dispersion 



relations. We have found that there is no significant difference in the massless tttx 
scattering rate when different dispersion relations are used. The rate of elastic 
scatterings is not affected by the speed of the pions. Assuming that pions are 
on the light-cone, to — k, the scattering contribution to the width of the pion of 



momentum k is studied at different temperatures. Figure |7.9| shows the results at 
T = 110 MeV and at T = 130 MeV. Harder momentum modes are damped more 
strongly. Damping of massless pions due to scattering turns out to be almost as 
significant as Landau damping. Pion modes thus can become broad. 

When T c is approached from above not only the heavy sigma fields but also 
the much lighter and even the massless pion fields are significantly damped. This 
fact should certainly be taken into account when considering formation and decay 
of DCCs. 
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Figure 7.6: Scattering contribution to the width of the pion as function of temper- 
ature in the model with explicitly broken chiral symmetry. 



115 



600 



500 
> 400 
S 300 
200 
100 



c mass 

C-7T 

a-a 
total 



20 40 60 80 100 120 140 160 180 200 

T [MeV] 



Figure 7.7: Scattering contribution to the width of the sigma meson as function of 
temperature in the model with exact chiral symmetry. 
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Figure 7.8: Massless pion-pion scattering rate versus momentum at T = 110 MeV 
(lower curves) and at T = 130 MeV (upper curves) with pion dispersion r elations 
J 1 = v 2 k 2 : v n = 1 (solid), v 2 = a (dashed) and vl = a + bk 2 (dotted). 
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Figure 7.9: Scattering contribution to the width of the pion as function of temper- 
ature in the model with explicitly broken chiral symmetry. 
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Chapter 8 
Conclusions 



In this dissertation different aspects of the non-equilibrium chiral/confinement 
phase transition as well as in-medium properties of the scalar and pseudoscalar 
mesons were investigated using quantum field theoretical descriptions at nonzero 
temperatures and densities. The SU(2)r x SU(2)l linear sigma model was used 
as an effective theory describing the chiral symmetry restoring phase transition of 
QCD. 

In Chapter the chiral phase transition at nonzero temperature and baryon 
chemical potential was studied at mean field level in the sigma model that includes 
quark degrees of freedom explicitly. Use of such a model was justified since high 
temperatures are dominated by quarks and low temperatures are dominated by 
mesons. The model was tuned to reproduce the correct properties in the vacuum. 
First the order of the phase transition was determined. For small bare quark 
masses we found a smooth crossover for nonzero temperature and zero chemical 
potential and a first order transition for zero temperature and nonzero chemical 
potential. The first order phase transition line in the (T, /i) plane ends in the 
expected critical point. The behavior of quark and meson masses with changing 
temperature and chemical potential was also examined. It has been found that the 
a mass is zero at the critical point. Adiabatic lines were computed. The results of 
our numerical calculations for the behavior of the adiabats shows a pattern which 
is in contrast with some earlier predictions on this subject [[Hp. In fact, the phase 
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transition found in this model turned out to be of the liquid-gas type. However, 
we showed that the strength of this transition depends strongly on the choice of 
coupling constants. 

The theory including quarks is well suited to temperatures above the critical 
temperature. In the rest of the work we kept our focus on the phase below the 
critical temperature. The model thus was slightly changed: quarks were dropped 
and fluctuations of meson fields were included beyond the mean field. 

In Chapter ^] the concept of disoriented chiral condensates (DCC) was reviewed 
in the context of their possible formation as a result of a non-equilibrium phase 
transition from a chiral symmetry restored phase to the symmetry broken phase. 
Coarse-grained equations of motion for the soft homogeneous and inhomogeneous 
chiral condensate fields were derived. The sigma and pion field configurations 
were coupled to a heat bath of hard momentum sigma mesons and pions. Based 
on linear response theory the fluctuations of the thermal bath, as a response to 
the presence of the non-thermal condensate, were analyzed. Multiple aspects of 
the response functions were emphasized and discussed in chapters that followed. 
Throughout this work the theory with spontaneously broken symmetry and the 
theory with explicit chiral symmetry bearking were examined separately. 

In Chapter |] the temperature-dependence of pion and sigma meson masses was 
studied. These mesons were dressed by the interactions with the thermal medium. 
Mathematically this was expressed via the response functions. A self-consistent 
evaluation of the meson masses including all one-loop order contributions was per- 
formed. The mass of the sigma meson is completely determined by the equilibrium 
condensate. Numerical results for the temperature-dependence of this show a first 
order chiral phase transition in the theory with spontaneously broken symmetry. 
However, we showed that changing the value of the coupling constant can render 
the transition second order. Special attention was given to the fulfillment of Gold- 
stone's Theorem when the chiral symmetry is spontaneously broken, according to 
which the pion is massless at all orders in perturbation theory and at all tempera- 
tures. In the model with explicit symmetry breaking a smooth crossover between 
the symmetry broken and symmetric phases has been found. 



120 



In Chapter |5|, the effect of the thermal medium on the velocity of long- wavelength 
pion modes was computed. Two different approaches were used to derive an ex- 
pression for the velocity of pions at finite temperatures. Calculating in coordinate 
space directly from the response functions and in momentum space via the ex- 
pansion of the pion self-energy about its small frequency and momentum yielded 
consistent results. Numerical analysis clearly showed a deviation of the velocity of 
massless pions from the speed of light at finite temperatures. Then the modifica- 
tion of the pion dispersion relation compared to the vacuum relation was discussed. 
The softening of the dispersion relation due to second-order contributions of the 
low momentum expansion was found. Going beyond this approximation and tak- 
ing into account fourth-order contributions resulted in a less drastic softening with 
respect to the vacuum relation. The effect was found to be increasingly impor- 
tant with temperature for a given momentum mode, and for harder momenta at 
a given temperature. The limit of validity of our model was determined by re- 
quiring causality not to be violated, that is, the group velocity should not become 
superluminous. 

Finally, in Chapters ^ and [7| the dissipation of the long-wavelength sigma and 
pion condensate fields was examined at one- and two-loop order, respectively. Dis- 
sipation of the condensate occurs because energy can be transferred between the 
condensate and the heat bath through the interactions of the different degrees of 
freedom. In Chapter |6] the effect at one-loop order was calculated by evaluating the 
response functions directly. Dissipation follows from first principles when deriving 
the equations of motion from the Lagrangian. The direct connection to the physi- 
cal processes responsible for dissipation was shown. In Chapter [7| the damping of 
the condensate fields due to elastic scatterings has been examined. These processes 
were neglected in the literature, based on the argument that they are higher order 
in coupling. Our results are contrary to previous beliefs according to which sigma 
mesons are quickly damped and, due to their small mass, pion oscillations persist 



in the small momentum region |[48|| . We found that at high temperatures not only 
is the damping of sigma mesons significant, but also that of the pions. This damp- 
ing is most accentuated in the phase transition region. The scattering rate of pions 
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becomes significant at high temperatures and it is even comparable in magnitude 
to the damping due to the one-loop order Landau damping. Dissipation of sigmas 
from the condensate due to their scattering was found to be much smaller than 
their dissipation due to their decay. The results of these two Chapters proved that 
dissipative processes up to two-loop order are significant when considering possible 
formation and decay of DCCs. 
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Appendix A 

Thermodynamics with Small 
Quark Masses 



In the chirally symmetric phase the constituent quark mass M q is small and it 
goes to zero in the chiral limit. Therefore, it is instructive to evaluate the thermo- 
dynamic potential Q for small M q . In this limit Q qq can be expanded in a power 
series in M q . The explicit expression is 



M 2 fd 2 Vt - 

n^T,KM q ) = WT,ti + —\—£ 



+ ... 



(A.l) 



M=0 



Here the first term VL (T, jj) = fl qq (T, /i, 0) can be easily calculated for arbitrary T 
and /i. The well-known result is 



2tt 2 



I^T 4 + ^TV + -V 
180 6 P 12 P 



(A.2) 



The quark number and entropy densities for massless fermions are obtained by 
differentiating Qq(T, /i) with respect to ji and T, respectively, 



n 



6tt 2 



n 2 T^ + ^) , 



(A.3) 



(A.4) 
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The second term in eq. ( |A.1| ) differs only by a factor of M q from the scalar density 
defined in eq. (|2.13 ). A straightforward calculation gives 



\ dM2 ) M = \MJm^o 




(A.5) 



This can be used to estimate the pion and sigma masses at large T and/or fi. 
Expressing M q in terms of mean meson fields, equation ( |2.7| ), and using the defi- 
nition of effective masses, equation ( [2.15| ), one arrives at the following asymptotic 



(M — > 0) expression for the pion and sigma masses 

Ml = Ml=g\{^- + ^-\ . (A.6) 



It shows that deep in the chiral-symmetric phase the pion and sigma masses are 
degenerate and large. At high temperatures T ^> \i and 

M^ = M a = gT , (A.7) 

where g ~ 3 in our calculations. Therefore, in this limit the contribution of 
pion and sigma excitations to the thermodynamical potential is suppressed by the 
Boltzmann factor er M °l T = ^T 9 . 
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Appendix B 

Imaginary Part of the Self-energy 



In the following we compute a very typical one-loop self-energy diagram, shown 
in figure B.l| , in which the bosons have a cubic coupling #00102. Here g is the 



coupling constant and 0, 0i and 02 have masses M, mi and m>2, respectively. 
Following general finite temperature Feynman rules |75| the contribution to the 



,p + k 




02, V 

Figure B.l: One- loop self-energy contribution to 
with cubic coupling to 0i and 02- 



self-energy of 0, with four-momentum k, in imaginary time (Euclidean space) is 
U(k) = U(u n ,k) = g 2 TY^ [ -^D l (uj m + uj n ,p + k)D 2 (uj m ,p) . (B.l) 



The free imaginary-time propagators have the form 

1 



D (u n ,p) 



w„ + p 2 + m 2 



(B.2) 
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and so the self-energy is 



U(u n ,k) = g 2 Tj2 J 



d 3 p 



{2^) 3 {u m + u n ) 2 + E 2 u 2 m + E 2 



2 - 



(B-3) 



where E\ = \J mf + (p + k) 2 and E± 



ri)-, 



p 2 are the on-shell energies of 



the internal lines, k is the momentum of the external boson and u n are the 
Matsubara frequencies that take only discrete values, u n = 2irnT, with n integer. 
It is worth noting that we keep the frequency u n and momentum k as independent 
variables, not requiring boson to be on mass-shell. There are standard ways of 
evaluating such Matsubara-sums. The interested reader is refered to ||7S|, [76] . My 
preference is replacing sums with contour integration. For a function f(p° = iuj n ) 
with no singularities along the imaginary p° axis the sum is given by 



T E /(P°) 



1 

2vri 
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dp \f(p°) + f(-p°) 
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(B.4) 



This replacement leads to the real time (Minkowski space) form of the self-energy: 



n(fc°,fc) = i 



d 3 p 
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k° — Ei — E 2 
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(B.5) 



1 1 

.k°-E 1 + E 2 ~ k° + Ex - E 2 

fi = f(Ei) = l/(exp (Ei/T) — 1) are the Bose-Einstein distribution functions of 
particle species i. One can extend this function to the whole complex plane by 
analytical continuation k° = iu n = u> + it. 

d 3 p 1 
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(2n) 3 AE 1 E 2 
+ (/ 2 
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uj — Ei — E 2 + ie uj + Ei + E 2 + ie 
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(B.6) 



,u — Ei + E 2 + ie u + Ei — E 2 + ie 

Since the physical interpretation of the various terms becomes more transparent 
in the imaginary part of the self-energy, in the following we derive the imaginary 
part of ( |B.6j ). Apply the identity known from complex analysis: 

1 „ 1 



k° - E±ie 



k°-E 



=F w5(k° - E) , 



(B.7) 
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where P means principle value. Then 



Imn(w, k) = -irg 2 J ^ f [(1 + h + -E*- E 2 ) - 5{u + E 1 + E 2 )) 

+(/ 2 - fx)(6(u -E 1 + E 2 ) - 5{u + E ± - E 2 ))] (B.8) 



B.l Direct Link to Physical Processes 

Based on relativistic kinetic theory one can write the rate for a process to occur 
as the square of the amplitude of the process weighted with the proper statistical 
factors, thermal distributions, and integrated over the phase space. Thus, the 
probability for a particle to propagate through aT / medium with energy 
u) < and momentum k will decrease with a rate V^ioj) and will increase with a 
rate T^u) given by 



T d (u, k) = 7^J dQ ab \M(<S>l...a l'...6)| 2 /i-/a(l + + h) (B.9) 

T t (uj, k) = ^J dQ ab \M(l'...b -> $l...a)| 2 /i'-/ 6 (l + /i)-(l + fa) (B.10) 
the phase-space 

dVl ab = dp 1 ...dp a dp 1 > ...dp b (2Tc) 4 5 4 k + ^Pi - X) Pj\ (B.ll) 

V *=i i=i' / 



dp = rf3 f ! (B.12) 
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